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Abstract 

These are pedagogical lecture notes on hydrodynamic fluctuations in normal relativistic 
fluids. The lectures discuss correlation functions of conserved densities in thermal equi- 
librium, interactions of the hydrodynamic modes, an effective action for viscous fluids, 
and the breakdown of the derivative expansion in hydrodynamics. 
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1 Basic hydrodynamics 



1.1 Introduction 

One can think about hydrodynamics as an effective long-distance description for a given 
classical or quantum many-body system at non-zero temperature. These lectures will be 
mostly about relativistic hydrodynamics, in other words, about hydrodynamics of fluids whose 
microscopic constituents are constrained by Lorentz symmetry, as happens in relativistic 
quantum field theories. Relativistic hydrodynamics is not limited of course to the description 
of the collective motion of particles that move at speeds close to the speed of actual light. All 
that matters is Lorentz symmetry: the hydrodynamic equations describing the quark-gluon 
plasma [1] will look identical to the hydrodynamic equations describing relativistic quantum 
critical regions in quantum magnets [2]. 

Besides its relevance to the real-world physical phenomena, hydrodynamics is an interest- 
ing subject from a purely theoretical physics point of view. It turns out that hydrodynamic 
fluctuations (both in normal fluids and in superfluids) can be thought of as gravitational 
fluctuations of black holes, and vice versa [31 HI- In fact, this connection between gravity 
and hydrodynamics works beyond the linear response, allowing one to derive how the full 
non-linear relativistic Navier-Stokes equations are encoded in Einstein's equations of General 
Relativity. This relation of the Einstein's equations to the equations of relativistic hydrody- 
namics is often referred to as the "fluid-gravity correspondence" [5]. 

What does this teach one beyond being a curious mathematical correspondence? One 
answer comes from the gauge-gravity duality [H [7j which states (among many other things) 
that the hydrodynamics of a black hole does in fact represent the hydrodynamics of a specific 
quantum system, whose Hamiltonian is explicitly known. In particular, one can easily evaluate 
transport coefficients in these quantum systems from black hole physics. This is an important 
advance because the quantum systems in question happen to be strongly interacting, and 
the conventional methods such as the diagrammatic expansion or Monte-Carlo simulations 
become inadequate to compute the transport coefficients. Because these quantum-mechanical 
models share important common features with the theory of strong interactions (QCD), the 
computations of transport coefficients from black hole physics becomes more than a purely 
academic exercise. 

An example of a fruitful interplay between hydrodynamics and black hole physics came 
from the realization that the kinematic viscosity in the gauge-gravity duality takes a universal 
model-independent value [HI [9], and that the Heisenberg uncertainty relation may prevent the 
existence of perfect fluids in nature [8] . The two main competitors for the most nearly-perfect 
fluid are the ultracold Fermi gases and the quark-gluon plasma, see [10] for a review. 

On the other hand, the relation between hydrodynamics and black holes has motivated 
new research into relativistic hydrodynamics itself, and its applications. One example is the 
classification of two-derivative terms (originally introduced in order to restore the causality of 
the standard one-derivative relativistic generalization of the Navier-Stokes equations) which 
only appeared recently [TH [T2] . Another example is the application of the relativistic hydro- 
dynamics to thermo-magnetic charge transport in high-temperature superconductors [13], and 
in graphene [14] . Another example is the realization that the standard relativistic hydrody- 
namics needs to be corrected in the presence of quantum anomalies |15j . Another example is 
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the modification of relativistic hydrodynamics due to parity violation [16]. Another example 
is the systematic study of dissipative relativistic superfluid hydrodynamics [TF] . All the above 
examples could have been understood many years ago without any black hole physics. But 
they weren't. Instead, the studies were heavily motivated by the recent experiments studying 
the quark-gluon plasma [18], by the gauge-gravity duality [191 120], and by the progress in 
understanding quantum critical phenomena [21]. 

These lectures will largely focus on hydrodynamic fluctuations, that is small, long-wave- 
length fluctuations near thermal equilibrium. These fluctuations are long-lived: their fre- 
quency w(k) vanishes as k — > 0, and they can propagate long distances, as one can see in the 
example of sound waves, whose dispersion relation is cu(k) = v s \k\ + 0(k 2 ). Such hydrody- 
namic modes are solutions to the hydrodynamic equations. At the same time, such modes are 
a danger to the very existence of hydrodynamics, as they violate the assumption of local equi- 
libration. Thermal excitations in these long-lived modes will contribute to both charge and 
momentum transport. For non-relativistic fluids, it has been known for a long time that this 
effect wipes off the very notion of classical hydrodynamics in 2+1 dimensions [22], and makes 
two-derivative corrections to classical hydrodynamics in 3+1 dimensions problematic [23] 1*1 




Relativistic hydrodynamics is no different in this respect: for small fluctuations near thermal 
equilibrium, there is not a great difference between relativistic and non-relativistic hydro- 
dynamics. Thus the classical equations of second-order relativistic hydrodynamics in 3+1 
dimensions do not describe the fluid correctly in the hydrodynamic limit [26] . 

These lecture notes have three main parts. The first part is about the equations of hy- 
drodynamics: we need to know what the equations are, before we solve them! The equations 
can be found for example in the book by Landau and Lifshitz [27], or in the book by Wein- 
berg [28J . The second part is about how the classical equations of hydrodynamics determine 
the two-point correlation functions of conserved densities, (energy density T 00 , momentum 
density T° l and charge density J°) in equilibrium. This will be done in the linear response 
theory, deriving the Kubo formulas for the shear viscosity, the bulk viscosity, and the charge 
conductivity. An excellent pedagogical reference on the linear response theory is the paper 
by Kadanoff and Martin [2S]. We will also compute the correlation functions by a somewhat 
different method, by studying the response of the hydrodynamic equations to the external 
gauge field, and to the external metric. The third part of the lecture notes is about the in- 
teractions of the hydrodynamic modes. In order to treat the interactions systematically, one 
may introduce microscopic random stresses and random currents, thus converting the classical 
hydrodynamic equations into stochastic equations. Using a standard method, hydrodynamics 
can then be recast as a real-time "quantum" field theory, where temperature plays the role of 
h, and where correlation functions follow from a generating functional defined by a functional 
integral with a local effective action. 

^dn systems with many degrees of freedom, these fluctuation effects can be suppressed if the limit N — > oo 
is taken before before the hydrodynamic limit k — > 0, uj — > 0. In other words, the large-iV limit does not 
commute with the hydrodynamic limit in 2+1 dimensions. In the context of the gauge-gravity duality, taking 
the limit N — > oo first is what allows the classical 2+1 dimensional hydrodynamics to emerge from 3+1 
dimensional gravity [531 12S] ■ 
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1.2 Non-relativistic hydrodynamics 

Let us start with non-relativistic hydrodynamics of normal fluids in three spatial dimen- 
sions. The discussion here will be a brief review; we will move onto relativistic fluids shortly 
which will be discussed in more detail. Hydrodynamic equations express conservation laws 
of whatever is conserved, supplemented by assumptions about the dissipative response, and 
thermodynamic information such as the equation of state. Let us first discuss ideal (i.e. 
non-dissipative) hydrodynamics of a simple non-relativistic fluid [27J. 

Ideal fluids 

In a non-relativistic system, particle number is conserved, in other words particles are not 
created or destroyed. This conservation of particle number is expressed in hydrodynamics as 
conservation of mass, by the continuity equation 

%> + ,%G<) = 0, (1.1) 

where p is local fluid density (mass per unit volume), and v is local fluid velocity. Another 
equation is the equation of motion of a fluid element, p(dvi/dt) = —d^p, where p is pressure, 
and (—dip) is the force per unit volume. In terms of the velocity field v(t, x) this gives the 
Euler equation, p(d t Vi + VkdkVi) = —dip. Using the continuity equation (11. ip . one can rewrite 
the Euler equation in the form of momentum conservation, 

d t (jnj i ) + d j n 4j = Q, (1.2) 

where Ily = pSij + pViVj. So far we have five unknown functions (p, v,p), and four equations. 
The equilibrium equation of state will provide a relation between p and p, but in general 
will introduce another unknown parameter such as temperature or entropy. This means that 
we need one more dynamical equation. The extra equation can be taken as the condition of 
entropy conservation, ds/dt = d t s + Vkdk5 = where s is entropy per unit mass. Equivalently, 
d t s + di(svi) = 0, where s is entropy per unit volume. Using the continuity equation, the 
Euler equation, and basic thermodynamic identities, the condition of entropy conservation 
can be rewritten as conservation of energy, 

d t (e+^)+a i {{w + ^) Vi )=0, (1.3) 

where e is (internal) energy per unit volume, and w = e + p is enthalpy per unit volume. We 
now have six unknown functions (p, v,p, e), and five dynamical equations (ll.ip -f FOj) . plus 
the equation of state which relates e,p and p. The above equations of ideal fluid dynamics 
express the conservation of mass, momentum, and energy. 

In writing down Eqs. (jl.ip - (jl.3p . we have implicitly assumed that the fluid is homoge- 
neous, for example there is only one species of particles. If there are several conserved "particle 
numbers" , the set of hydrodynamic equations needs to be expanded, reflecting the fact that 
there are more conserved quantities. See for example |27], §58. 
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Figure 1: A stationary flow of an ideal fluid with a velocity gradient. Thin arrows repre- 
sent particles which can transfer x-momentum in the y-direction, eventually leading to the 
equilibration of the inhomogeneous velocity profile. 

Dissipative fluids 

Equations (ll.ip — (jl.3p of ideal hydrodynamics ignore dissipation. For example, sound waves 
derived from ideal hydrodynamics would propagate indefinitely without damping. Physically, 
dissipation means that the fluid can move locally not only because a fluid element is pushed 
around by other fluid elements, but also because the microscopic movement of the constituents 
of the fluid can even out various inhomogeneities in velocity, temperature etc. For example, 
consider the configuration with constant p, p and e, and with 



This describes homogeneous fluid which moves uniformly in the x direction, but whose velocity 
depends on y, as shown in Fig. CD For such configuration, = 0, as well as v^d^Vi = 0, and 
therefore the flow ( 11. 4ft is a perfectly good stationary solution to the ideal hydro equations 
( ll.ll) -( fl~3l) . However, this flow is unphysical because the particles which comprise the fluid 
can move in the y direction, and this will lead to momentum transfer between different layers 
in the fluid. As a result, the velocity distribution v x (y) will tend to become more uniform, 
and the configuration (11.41) will not remain stationary. It is precisely these kind of processes 
that are taken into account by the dissipative terms in the hydrodynamic equations. The 
idea of dissipative hydrodynamics is to modify the equations of the ideal hydrodynamics by 
adding extra terms which are proportional to spatial gradients of the hydrodynamic variables. 
Moreover, one assumes that the gradients are small, so that a "gradient expansion" can be 
constructed. In hydrodynamics one always assumes that the quantities of interest vary on 
scales which are large compared to the microscopic length scale Z m f p , so that \dip\ <C |p// m f P |- 
We will retain only those terms which are linear in spatial gradients. In principle, one can 
write down dissipative terms which are quadratic in spatial gradients, then cubic and so on. 
One can think about hydrodynamics as of an effective theory which is valid at long distance 
and time scales, similar to the low-energy chiral theory in QCD. Higher-order terms are 
suppressed by the powers of the cutoff which in hydrodynamics is Z m f p . 



v x = v x (y) v y = v z = 0. 



(1.4) 
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The equations of dissipative hydrodynamics can be written in the form [27] 



d t p + di(pvi) = , (1.5a) 

dtifWi) + djIUj = 0, (1.5b) 

^(e + ^)+^ = 0, (1.5c) 

where the stress tensor is Uij = p5ij + pViVj — with 



Ey = rj (diVj + djVi - -6ij d k v k ) + CSijd k v k 



and the energy current is 



j- = ( w+^- ) Vi - YjijVj - ndiT . 



The coefficient rj is called shear viscosity, £ is called bulk viscosity, and k is called thermal 
conductivity!^ The transport coefficients rj, (, and k should be thought of as input parameters 
for hydrodynamics: like in any effective theory, they can not be evaluated in hydrodynamics 
itself, but rather need to be computed from the underlying short-distance physics, for example 
using the Boltzmann equation |31] . In a near-ideal gas, both r\ and k are proportional to the 
mean-free path. 

In writing down the equations (11.51) of dissipative hydrodynamics, a particular out-of- 
equilibrium definition of the hydrodynamic variables has been chosen. A different choice 
would lead to the hydrodynamic equations which look different from ( 11. 5ft . We will discuss 
this issue when talking about relativistic hydrodynamics in Section 11.31 

One final comment concerns the terminology. The words "perfect fluid" , "ideal fluid" , or 
"ideal hydrodynamics" usually refer to non-dissipative hydrodynamics, i.e. the form of the 
hydrodynamic equations with rj, (, and k set to zero. However, ideal hydrodynamics is not 
the same as hydrodynamics of the ideal gas. In fact, ideal gas (i.e. gas without interactions 
between the particles) has no hydrodynamics at all because the thermal equilibrium state is 
unstable when there are no interactions. This pathology of the ideal gas can be cured by 
introducing small interactions between the particles, with interaction strength determined by 
a small parameter A. The equation of state in such a system is the same as the equation of state 
of the ideal gas, up to small corrections of order A. The mean free path, on the other hand, 
will be very large (inversely proportional to A), and the length scale at which hydrodynamics 
is applicable must be even larger. The shear viscosity and the thermal conductivity in this 
slightly non-ideal gas are large, because rj and k are directly proportional to the mean-free 
path. By taking the limit A — > one can see that the ideal gas has infinite shear viscosity 
and thermal conductivity, unlike the ideal fluid which has (by definition) zero shear viscosity 
and thermal conductivity. 



2 Often rj is called just "viscosity", and £ is called "second viscosity" or "volume viscosity". The equation of 
momentum conservation (jl.5bl) is called the Navier-Stokes equation. For a historical account of the formulation 
of the Navier-Stokes equation, see [30] . 
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Hydrodynamic modes 

Now that we have the equations of hydrodynamics, we can study some simple solutions. A 
fluid at rest, with constant p=p, p=p, e=e, and v=0 is obviously a solution. Now let us look 
at small fluctuations around this solution: p=p + Sp, e=e + Se, p=p + Sp, and small v. "Small 
fluctuations" means that Sp <C p etc, but we need to say what does it mean to have small v. 
We will see that the relevant condition is that v must be much smaller than the speed of 
sound. 

Let us start with ideal hydrodynamics. The linearized equations (ll.ip - (11.3j) become 

d t 5p + pdiVi = 0, (1.6a) 
p d t Vi + di5p = , (1.6b) 
d t 8e + w diVi = , (1.6c) 

where w = e+p. There are three unknown functions besides v, and we need to choose 
the variables: given the equation of state, we can work with Se, 5p, or with Sp, ST, or with 
some other variables. Let us choose Se and Sp as the hydrodynamic variables, then we need 
an equation of state of the form p = p(e,p). In practice these variables may not be the 
most convenient ones, but they are natural from the effective field theory point of view 
because they represent the densities of conserved charges. From the equation of state, we have 
d{Sp = (dp/de) p diSe + (dp/dp) e diSp, where the thermodynamic derivatives are determined 
by the equilibrium equation of state. 

For fluctuations around an infinite homogeneous background, we can take all variables 
proportional to e* kx , and decompose the velocity as v = + v ± , where the parallel and 
perpendicular are with respect to the wave vector k. In particular, v n = k(k-i;)/k 2 . Then 
it is easy to show using Eqs. (11 .6p that the transverse velocity component becomes time 
independent and decouples. On the other hand, the longitudinal component satisfies the 
wave equation 

af« ll (k,t) + k\ 2 « ll (k,t) = o, 

where v 2 s = (dp/dp) e + (dp/de) p w / p. This describes a hydrodynamic mode which propagates 
with the velocity v s . As is evident from Eqs. (11.61) . both mass density and energy density are 
part of the same mode, and describe density fluctuations which propagate with the velocity v s . 
This hydrodynamic mode is of course the sound wave, and v s is the speed of sound. The sound 
waves are longitudinal because v ± is not a part of the propagating mode. Also, Eqs. (I1.6P say 
that \Sp/p\ = |V||/f s | an d \Se/w\ = \v }l /v s \, so the magnitude of velocity fluctuations must be 
much smaller than the velocity of sound in order for the linearized approximation to be valid. 
The eigenfrequencies of the system (11.61) are w(k) = ±|k|t> s , plus two uj = solutions, one for 
v ± , and one for (Sp, v^). Once dissipation is included, the u = solutions will become the 
transverse non-propagating shear mode and the longitudinal non-propagating heat mode. 

Let us now bring the speed of sound into a more familiar form. To evaluate the thermody- 
namic derivatives, we need to choose an equilibrium ensemble. Let us work in the canonical 
ensemble with a fixed total number of particles N. Mass density is p = mN/V, where m is 
the mass of one particle. At fixed N, we have TdS = dE + pdV ', as well as dV/V = —dp/ p. 
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Also, the change in energy is dE/V = d(eV)/V = de — edp/p, and thus 



TdS 



de 



w dp 



V p 

This gives the variation of pressure in terms of dS and dp: 



TdS 



+ 



dp\ w /dp 



dp . 



The coefficient in the square brackets is precisely t> 2 , and therefore 

dp\ 
d PJs,N 

which is the familiar textbook expression. For an ideal gas, p = pT/m, and also p oc p 7 , 
where 7 is the adiabatic index, about 1.4 for air. So the speed of sound in a non-relativistic 
ideal gas is 

m ' 

proportional to the square root of temperature. 

The same analysis of small fluctuations around a static thermal equilibrium state can be 
repeated for the dissipative hydrodynamics. The dispersion relations w(k) for the hydrody- 
namic modes will get imaginary parts proportional to k 2 (at small |k|) which are determined 
by the dissipative coefficients r], ( and k. As a result of the 0(k 2 ) dissipation, low-frequency 
sound waves propagate farther than high-frequency sound waves. 



1.3 Relativistic hydrodynamics 

Now let us move on to relativistic hydrodynamics. As it is the main subject of these lectures, 
we will discuss it in some more detail. For simplicity, we will only discuss normal fluids (no su- 
perfTuids), and fluids without dynamical electromagnetic fields (no magneto- hydrodynamics). 
The speed of light is set to one. 



Hydrodynamic variables 

Just like in the non-relativistic case, hydrodynamic equations express conservation laws of 
whatever is conserved. According to the Noether theorem, conservation laws are related to 
continuous symmetries of the fundamental microscopic theory, which imply the existence of 
conserved currents. In relativistic systems, the spacetime symmetries are given by transla- 
tions, rotations, and boosts. In addition, there may be other "internal" symmetries, such as 
the U(l) baryon number symmetry. 

The conserved current corresponding to the spacetime translation symmetry is the energy- 
momentum tensor T^ v . We will define the energy- momentum tensor by the variation of the 
action with respect to the external metric [28], in which case it is symmetric, T^ v = T u ^, and 
the conservation law is 

8^ = 0. (1.7a) 
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The conserved currents corresponding to rotations and boosts are Ai tluX = x fl T uX — x u T^ lX , 
which are identically conserved owing to the symmetry of the energy- momentum tensor, 
d\A4 flvX = 0. In other words, given our definition of the energy-momentum tensor, there is 
no extra conservation equation beyond (11.7aj) expressing the symmetry under rotations and 
boosts. Further, we will allow for a possibility of a conserved current J M corresponding to a 
global U(l) symmetryjf) 

<V M = °- ( L7b ) 
In d spatial dimensions, there are d+1 equations (11.7aj) . and one equation (11.7bl) . while there 
are (d+l)(d+2)/2 components of T^, and d+1 components of J M , in other words there are 
more unknowns than there are equations. The simplifying assumption of hydrodynamics is 
that T^ v and J M can be expressed as functions of d+2 fields: a local temperature T(x), a local 
fluid velocity v(x), and a local chemical potential fi(x). This makes the number of unknowns 
equal the number of equations. 

The choice of hydrodynamic variables can be understood as follows. The equilibrium 
state is characterized by the density operator g, which is proportional to the exponential of 
the conserved charges [32] . 

g = ^e^ p ^ N . (1.8) 

In this expression Z = tr e^* +7 is the partition function, is the momentum operator, 
and N is the operator of the conserved charge. The density operator is characterized by a 
timelike vector /3 M and a scalar 7, which parametrize the manifold of thermal equilibrium 
states. We will write these parameters as (3^ = ftu^, and 7 = where the vector is 
normalized such that u 2 = —1, and represents the velocity of the fluid, — 1/T is the inverse 
temperature, and fi is the chemical potential. In components, = (1— -y 2 ) _1 / 2 (l, v), where v 
is the spatial velocity. Note that the equilibrium state breaks Lorentz symmetry (specifically, 
boost invariance) due to the presence of the preferred timelike vector 

Hydrodynamic is concerned with states which deviate slightly from the thermal equilib- 
rium states (II. 8p specified by constant /3 M and 7. Hence, one chooses to describe these slightly 
non-equilibrium states in terms of slowly varying functions fl^x) and 7(2), or equivalently 
in terms of slowly varying functions u^(x), T(x), and n(x). As the microscopic theory is 
Lorentz invariant, and it is only the state (I1.8P that breaks the Lorentz symmetry, the hydro- 
dynamic equations formulated in terms of (3^(x) and 7(2;) must be Lorentz covariant, and the 
effective action for hydrodynamics must be Lorentz invariant. We will return to the question 
of effective action in Section HI In the following, we will focus on equations of relativistic 
hydrodynamics, and will use u^, T, and \i as our hydrodynamic variables, as is conventional. 
One could have used the variables (3^ and 7 instead. 

Constitutive relations 

For any timelike vector u M , the energy-momentum tensor and the current can be decomposed 
into components which are transverse and longitudinal with respect to w M , using the projector 



3 For simplicity, we will assume a parity-symmetric microscopic theory, and take to be a vector current 
(charge density is parity-even). If parity is not a symmetry of the microscopic theory, or if J M is an axial 
current (charge density is parity-odd), there will be extra terms in the constitutive relations in addition to 
those we are going to discuss here [17l [16] . 



9 



A^ u = rf u + u' 1 u u , where rj^v = diag (— 1, 1, 1) is the flat-space metric. We write the 
decomposition as [33] 

= Su^u" + VA^ + {<fu v + gV) + iT , (1.9a) 

J* = MvP + f, (1.9b) 

where the coefficients £, V, and M are scalars, the vectors and j M are transverse, i.e. 

= = 0, and the tensor t^ u is transverse, symmetric, and traceless. Explicitly, the 
coefficients are: 



1 

1 



£ = u„u u T^, P = -A ia ,T"", M = - Utl J\ (1.10a) 



g„ = —A m upT ap , j„ = A^ V J U , (1.10b) 
V = \ UVaA^ + A ra A M/3 - ^A^aJ) T afi . (1.10c) 

Equations (II. 9ft are just identities which hold locally in spacetime for any symmetric T^ v [x) 
and any J^{x), given a chosen vector u^(x). The assumption of hydrodynamics enters in 
expressing the coefficients £, V, q^, , A", and in terms of the hydrodynamic variables 
Up, T, and \i. For example, the coefficients £, V, and A" will be functions of scalars T, /i, c^-?/ 1 , 
u^d^T, d 2 fi, (d^u u ) 2 etc. Similarly, the coefficients and q^ will be functions of transverse 
vectors A fUl d v T ', A fiu d u fi etc. The expressions for T^ u and J M in terms of the hydrodynamic 
variables T, tt M , and \i are called the constitutive relations. 

We will write down the constitutive relations in the derivative expansion, which is the ex- 
pansion in powers of derivatives of the hydrodynamic variables. Because the deviations from 
thermal equilibrium are assumed small, one expects that one- derivative terms will be smaller 
than zero-derivative terms, two-derivative terms will be smaller than one-derivative terms and 
so on. This is analogous to the derivative expansion in effective field theory, while the role of 
the cutoff is played by the microscopic scale l m f p . At length scales shorter than i m f p , hydro- 
dynamics stops being a sensible description of matter. Ideal (non-dissipative) hydrodynamics 
corresponds to taking into account only non-derivative terms of the hydro variables T, u^, 
and fi in the expansion of T^ v and J M . Conventional dissipative hydrodynamics corresponds 
to taking into account both non-derivative and one-derivative terms of the hydro variables T, 
Up, and // in the expansion of T^ u and J M . 

Zeroth-order hydrodynamics 

Let us start with ideal hydrodynamics, i.e. no derivatives of the hydro variables. The trans- 
verse q^, t^", and can only be built out of derivatives of the hydro variables, hence there 
are no q^, t^, and j M in the ideal hydrodynamics. On the other hand, the coefficients £, V, 
and M can be functions of T and fi. These coefficients have simple physical interpretation: 
in static equilibrium, T^ u = diag(e,p, ...,p), where e is the equilibrium energy density and p 
is the equilibrium pressure, while J M = (n, 0), where n is the equilibrium charge density. For 
a fluid which moves with constant velocity u^, the energy-momentum tensor and the current 
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are obtained by performing the corresponding Lorentz transformation, and one finds 



r"" = e«V+pA^, (1.11) 
J» = nu>*. (1.12) 

Ideal hydrodynamics corresponds to adopting this equilibrium form of the energy-momentum 
tensor and the current, and promoting e, p, u^, and n to slowly varying fields. The constitutive 
relations in ideal relativistic hydrodynamics thus have the form (II. lip . (11.12[) . in other words 
we can identify £(x) = e(x) with local energy density, V{x) = p(x) with local pressure, 
J\f(x) = n(x) with local charge density, and u^{x) with local fluid velocity. The equilibrium 
equation of state provides one with p(T, p), from which the energy density e, entropy density s, 
and the charge density n can be found as s = dp/dT, n = dp/ dp,, e = —p + Ts + pn. 
The longitudinal component of Eq. (I1.7ap . u v d^T^ u = 0, and current conservation give 

d tt {wu ft )=u' i d tl p, (1.13a) 
<9 M (ra^) = 0, (1.13b) 

where w = e + p is the density of enthalpy. Using the thermodynamic relation w = Ts + pn 
together with dp = s dT + ndp,, the conservation equations in 1 11. 13ft combine to give 

d fl (sun = 0. 

This is interpreted as the conservation of the "entropy current", namely that (locally) the 
entropy does not increase in ideal (non-dissipative) hydro dynamics 

First-order hydrodynamics: frame choice 

Let us now proceed to first-order hydrodynamics, i.e. taking into account terms in the con- 
stitutive relations with up to one derivative of the hydrodynamic variables. This brings up 
a subtlety which was not apparent in the ideal hydrodynamics. Namely, the notion of local 
temperature, local chemical potential, and local fluid velocity are not uniquely defined out 
of equilibrium. One may define many local temperature fields T(x) which will differ from 
each other by gradients of the hydro variables, and will approach the same constant value 
in equilibrium, when the gradients tend to zero. The same applies to the chemical potential 
p(x) and to the velocity u^(x). In particular, this means that coefficients E, V, and M in the 
decomposition (11 .9p are 

S = e(T,p) + fs{dT,dp,du), (1.14a) 
V = p(T,p) + f r (dT,dp,du), (1.14b) 
J\f = n(T,p) + f N (dT,dp,du), (1.14c) 

where e, p, and n are determined by the equation of state in equilibrium, and the form of 
the out-of-equilibrium gradient corrections fg, f-p, fx depends on the definition of the local 
temperature, local chemical potential, and local velocity. In hydrodynamics, redefinitions of 



4 Even in zeroth-order hydrodynamics, entropy can increase in presence of shock waves which are disconti- 
nuities of the hydrodynamic flow |27j . 
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the fields T(x), fi(x), and u^{x) are often referred to as a choice of "frame", and we adopt 
this terminology here. When thinking about frame choices, it is important to keep in mind 
that the parameters T(x), fi(x), and u^ix) have no first-principles microscopic definition out 
of equilibrium. For example, there is no "temperature operator" whose expectation value 
in a given non-equilibrium state would give a local temperature T(x). Rather, T(x), 
and u^(x) are merely auxiliary parameters used to parametrize T^ u (x) and J M (x), which do 
have microscopic definitions. The auxiliary parameters may be redefined at will, as long as 
the energy-momentum tensor and the current remain unchanged under such redefinitions, see 
e.g. [17] for a recent discussion. 

Consider a redefinition (or "frame" transformation) 

T(x) ->■ T\x) = T(x) + 5T(x) , 
fi(x) — » n\x) = n(x) + 8fi(x) , 



u 



where ST, 5fi and 5m m are first order in derivatives. Note that 5u^ must be transverse, 
u^Su^ = 0, due to the normalization condition u 2 = —1. Using the definitions (ll.lOj) . and 
remembering that q^, and t^ v are first order in derivatives, while T^ v and J M remain 
invariant, one finds to first order 

5£ = 0, 5V = 0, 6N=0, (1.15a) 
% = -(S+V)6u tl , 5 JfM = -MSu^ , (1.15b) 
5V = 0. (1.15c) 

One can choose 5u^ such that = 0; this definition of the local fluid velocity which is adopted 
by Eckart [23], is often referred to as "Eckart frame", and implies no charge flow in the local 
rest frame of the fluid. Or one can choose 5u^ such that = 0; this definition of the local 
fluid velocity which is adopted in the book of Landau and Lifshitz [27J, is often referred to 
as "Landau frame" , and implies no energy flow in the local rest frame of the fluid. Equation 
(jTT5a|> says that e(T,fi) + f £ (dT,dn,du) = e(T',/i') + f £ {dT' ,dii' ,du'), and similarly for V 
and jV, which implies that 

This means that we can use the redefinitions of T and \x to set two of the three functions 
f £ , f v , f^- to zero. It is conventional to set f' £ and fj^ to zero, which is to say that one can 
choose an off-equilibrium definition of T and /i such that £ = e and M = n. 

Other frames choices are possible: a generic frame would have both q^ and non-zero, 
as well as f £ , f-p, fx (defined by Eq. (11.140 ) all non-zero. As is evident from Eq. (11.15bj) . the 
combination 

71 

I* = f ff 

e + p 
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is frame-invariant. Using the above transformations of fg, f-p, and f^, it is easy to check that 

f = fr- (|) a - (I)/* 

is frame-invariant as well. The dissipative terms may be shifted around in the constitutive 
relations by frame transformations in a way that leaves B 1 and / invariant. 

First-order hydrodynamics: constitutive relations 

We start with the Landau frame, where the velocity is chosen so that g M = 0, and T and 
H are chosen so that E — e and M = n. The remaining coefficients V, t^", and in Eqs.f ll.9j ) 
need to be expressed in terms of the hydrodynamic variables, providing constitutive relations 
in first-order hydrodynamics. 

There are three one- derivative scalars in addition to T and /i, which can be built out of the 
hydro variables: u x d\T, u x d\fi, and d\u x . There are also three transverse vectors, A^ v d v T , 
A^c^/i, and A^u^, where u u = u x d\u v . Finally, there is one transverse traceless symmetric 
tensor, 

= A m A vfi (d a Uf, + dpU a - 2 -T) a pdX 

Let us look at the scalars first. We would write the derivative expansion as 

V = p + c lU x d x T + c 2 u x d X fi + c 3 d x u x + 0(d 2 ) , 

where p is thermodynamic pressure in the local rest frame of the fluid, 01,2,3 are some co- 
efficients, and 0(d 2 ) denote the terms which contain either second derivatives of the hydro 
variables, or are quadratic in the first derivatives. We can simplify our life by making use 
of the known zeroth-order hydro equations: indeed, two scalar equations u^d v T^ u = and 
d^J^ = provide two relations among the three one- derivative scalars. We can use the 
zeroth-order hydro equations to eliminate two of the scalars in the above expansion for V, 
and the error we are making by using the zeroth-order equations is only going to enter at 
0(d 2 ), which we are neglecting anyway. It is conventional to eliminate u x d\T and u x d\fi, and 
keep d\u x as the only independent scalar at first order. Thus the expansion of V becomes, 

V = p-(d x u x + 0(d 2 ), 

where ( is a coefficient which needs to be determined from the microscopic theory. The 
coefficient ( is the bulk viscosity. 

Similarly, we would write as a combination of the three transverse one-derivative vectors. 
Again, we can use the zeroth order hydro equations: there is one transverse vector equation 
Ax^d^T^ = 0, and therefore only two of the three transverse vectors are independent. We 
choose to eliminate A^ v u v , and therefore the expansion of j M becomes 

f = -oTA^dvQjt/T) + XtA^T + 0{d 2 ) , 

where a and Xt are coefficients which need to be determined from the microscopic theory. 
The coefficient a is the charge conductivity; we will see later that the coefficient Xt must be 
zero. 
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Finally, a^ v is the only transverse traceless symmetric tensor, hence 

= - V a» v + 0(<9 2 ), 

where the coefficient rj is the shear viscosity. Thus the constitutive relations in first-order 
relativistic hydrodynamics in the Landau frame can be taken as: 

T» v = eu"u v + pA"" - r] A>* a A vfi \d a up + d p u a - ^Va^uA - (A^d x u x + 0(d 2 ) , (1.16a) 

= nu^ - aTA^d v (/i/T) + X tA^8 u T + 0{d 2 ) . (1.16b) 

We see that Lorentz covariance restricts the constitutive relations up to four transport coeffi- 
cients rj, (, a, and Xt- We will see later that rj, £, and a must be non-negative, while Xt has 
to vanish (even though a non-zero value of Xt is allowed by Lorentz symmetry). The equation 
of state provides one with p(T, /x), from which one can find e(T, /x) and n(T, /x). The transport 
coefficients 77, £, and a depend on T and /x in a way that is determined by the underlying 
microscopic theory. 

In a generic unspecified frame, the constitutive relations for frame-invariant quantities in 
first-order hydrodynamics take the form 

iT = -rj , f = -C d x u x , r = -aTA^d v (ii/T) + X T^d u T . 

For example, one could choose a frame in which the bulk viscosity appears as a non-equilibrium 
correction to charge density, and the charge conductivity appears in the constitutive relations 
for both the energy- momentum tensor and the current. In the Eckart frame, the constitutive 
relations become 

= eu^ + pA^ + {q^u v +q u u^) -ija 1 " ' - (A^dxu* + 0{d 2 ) , (1.17a) 
r = nU » + 0{d 2 ) , (1.17b) 

where = {aTA^ u d v {ix/T) — XTA^ v d v T){e+p)/n. The expression for g M can be rewritten in 
an equivalent form by eliminating A^c^/x in favor of A^ v d v T and A^ u u V) using the equations 
of ideal hydrodynamics: 

g M = - K A» v {Tu v + d v T) - € -^- Xr A^d v T + O{0 2 ) , 

n 

where k = o~(e+p) 2 /(n 2 T) is the heat conductivity. It is this Eckart-frame form of the 
relativistic hydrodynamic equations with % T = that becomes Eq. (11. 5p in the non-relativistic 
limit. 

The coefficients rj, (, a, and Xt can be viewed as parameters in the effective theory 
that need to be matched to the microscopic theory. This matching can be done in the linear 
response theory, as we will see in the next section. The linear response theory will give explicit 
expressions for 77, (, and a in terms of correlations functions of the energy-momentum tensor 
and the current. As the latter are independent of the "frame", this provides an alternative 
way to see that the transport coefficients rj, (, and a are frame invariant: choosing a different 
frame may change the place where they appear in the constitutive relations, but not their 
value. 
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As this example of first-order hydrodynamics illustrates, Lorentz covariance alone is not 
sufficient to enumerate the transport coefficients: the number of non-vanishing transport 
coefficients does not coincide with the number of "independent tensor structures" that one can 
write down consistent with Lorentz symmetry, up to a given order in the derivative expansion. 
Other ingredients are needed (such as the local form of the second law of thermodynamics, or 
the linear response theory) in order to find the number of independent transport coefficients. 
The same happens in second-order hydrodynamics [3U ESS 136] - 

First-order hydrodynamics: entropy current 

Before moving on to the linear response theory, let us briefly discuss the entropy current 
in first-order relativistic hydrodynamics. In thermal equilibrium state with constant the 
entropy current is = su^. We saw that d^su^ 1 ) = in zeroth order (ideal) hydrodynamics. 
One now assumes that in first-order hydrodynamics there exists a current 

5 M = su 11 + (gradient corrections) , 

which satisfies d^S^ #J when the equations of hydrodynamics are satisfied. The gradient 
corrections vanish in equilibrium, and are built out of derivatives of T, p, and u^. The entropy 
current defined this way is in general not unique, however, demanding that d^S^^O, one can 
obtain the constitutive relations together with constraints on the transport coefficients [271137]. 
An expression for the entropy current may be written down by using a covariant version of 
the thermodynamic relation Ts = p + e — pn as [32] 

TS^ = pu» - T» u u u - pJ". 

Using the decomposition of the energy-momentum tensor and the current fll.9f> . we find 



5" 



u^ + ^-^f, (1.18) 



where £ , Af, q^ and are defined by fll.9p and (jl.lOp . Using frame transformations fll . 15j) 
together with Sp = p(T', p') — p(T, p) = s5T + n8p, it is straightforward to check that the 
entropy current fll . 18j) is frame invariant. The expression for the entropy current simplifies in 
a frame with £ = e, M = n, and in addition one can use either the Eckart frame condition 
= 0, or the Landau frame condition = 0. The positivity of entropy production d^S^ 1 ^ 
will dictate that 

00, C^O, (O0, Xt = 0. (1.19) 

The canonical form f ll . 18j) of the entropy current and the constraints fll. 19jl arise by demanding 
the positivity of entropy production in curved spacetime, V^S^ = 0, see Ref. [TTJ. We will 
see that the same conditions (I1.19P can be obtained in the linear response theory, without 
reference to an entropy current. 



2 Hydrodynamic equations and correlation functions 

If we wish to study small hydrodynamic fluctuations about a static thermal equilibrium state, 
we can take the full non-linear relativistic hydrodynamic equations fll.7[) with constitutive 



15 



relations (|1.16p . and linearize them around the solution u M = (1,0), T = const, fi = const. 
For most of this section, rather than working with T, u\ and /i, we will work with their 
conjugate variables, which are the energy density, momentum density, and charge density. 
The reason is that the latter have a microscopic definition given by the operators T 0/i (x) 
and J°(x). For a pedagogical discussion of the linear response theory in non-relativistic 
hydrodynamics, see the beautiful paper by Kadanoff and Martin [29] . 



2.1 Simple diffusion 

Before moving on to the full linearized hydrodynamics, we start with a simple example of the 
diffusion equation for a conserved density n(t,x), 

d t n-DV 2 n = 0, (2.1) 

where V 2 = d % di is the spatial Laplacian. This equation of course describes the diffusion of 
n(t, x) with a diffusion coefficient D, and gives a simple example of relaxation with u(k) — > 
as k— >0. To proceed, we define the spatial Fourier transform by 

" (t ' x) = /(2^ elk ' X " (t ' k) - 

We can now solve the equation, 

n(t,k) = e- z?k2 'n (k) , (2.2) 

where rio(k) = n(t=0, k). What this equation tells us is how the density relaxes at t > given 
an initial disturbance no(k). We will be working later with functions of frequency, rather than 
time, and so we would like to Fourier transform (12. 2\\ in time as well. However, n(t, k) in 
Eq. (12. 2p is only defined for t ^ 0, so instead of the Fourier transform, we will use the Laplace 
transform in time, defined as 



oo 

izt„ 



Jo 

where z must have a positive imaginary part for convergence. Applying the Laplace transform 
to the diffusion equation (12. II) . we find the solution as 

n(*,k) UM 



-iz + L>k 2 



We can also write this in terms of the source fi(t, x) for the charge density, which is the 
chemical potential. For small fluctuations, we have n(t, x) = x/i(t, x), where x = (dn/d/j,)^ 
is the static susceptibility. Note that x is a static thermodynamic quantity, not a dynamic 
response quantity. We now write the solution to the diffusion equation in terms of the initial 
value of the source as 

<-) 

We now would like to connect the solution to the diffusion equation to correlation functions 
in our system. To do so, we imagine the following process. We turn on a source for the density 
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of the form fi(t, x) = e £t fi(x)8(—t). In other words, the source /x(x) is adiabatically turned 
on at t = — oo, and then switched off at t = 0. We then watch the system evolve at time 
t > 0. The master formula for the whole of linear response theory comes from first-order 
time-dependent perturbation theory in Quantum Mechanics. For a system with a time- 
independent Hamiltonian H and a Heisenberg-picture operator A(t, x) defined with respect 
to H, we imagine adding a time- dependent contribution 5H(t) to the Hamiltonian. Then the 
change in the expectation value of A(t, x) is 

6(A(t, x)> = -i /" dt' ([A(t, x), 5tf(t')]> , (2-4) 

J — oo 

to first order in In our case 5H is due to the external source /x(t, x), and has the form 

5H = - d d x fi(t, x) n(t, x) 



[think grand canonical density operator, if — >• H—fiQ]. Now, from Eq. (12.41) . the change in 
the density is 



(n(t,x)) = i I dt' e £t '6{-t') I d d x' n(x!) ([n(t,x),n(t',x')]) 

J — oo J 

= i I dt' e £t 'd(t-t') I d d x' /i(x') ([n(t, x), n(t', x')]) . 

This can be written more compactly if we introduce the retarded function 

G£ n (t-t',x-x') = -i6{t-t') ([n(t,x),n(t',x')]) . 
Then the change in density induced by the external source fi(t,x) becomes 

(n(t, x)) = - f dt' e £t ' I d d x' /i(x') G* n {t-t', x-x') . (2.5) 

J — oo •/ 

Our goal now is to take this equation for t > 0, find the Laplace transform (n(z, k)), and 
compare with the prediction (12. 3ft of the diffusion equation. This will give us the retarded 
function. 

We first Fourier transform Eq. ( 12.51) in space, which gives 



(n{t, k)> = - f dt' e £t 'fi(k) G* n (t-t', k) . 

J — oo 



(2.6) 



Note that /x(k) is the Fourier transform of /i(x), in other words it is the Fourier transform of 
the external source x) at t=0, hence we will write /i(k) = /io(k), to use the notation of 
Eq. (12.31) . Next, we Fourier transform the retarded function in time, 



G R (t-t',k)= r — G R (u,k)e" 
J-oo2vr 



■iu>(t-t') 
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Note that G R (t, k) is only non-zero for t > 0, hence G R (u, k) is an analytic function in 
the upper half-plane of complex u. We can then analytically continue G R (u, k) to lower 
half-plane. Now the density induced by the external source becomes 



(n(t,k)) = -Mk) J^G n >,k); 



<'<Lc . ) e~ iuJt 

ioj + s 

We multiply both sides by e lzt (with lm.z > 0), and integrate over t from to oo, which gives 



2tt (ico + e) (i(u—z) + e) ' 

To do the integral, we close the contour in the upper-half plane where G R is analytic. There 
are two poles inside the contour, at uj = is, and co = z + is, thus 

(n(z, k)> = -/i (k) G '"" ( ' :,k) ~ G "" ( ' :=0,k) , (2.7) 

iz 

where the argument of G R is understood to be slightly above the real axis. Now comparing 
with equation f)2.3p . we find: 

G R n (z,k)-G R n (z=0,k)- 



-iz + Dk 2 

The only missing piece is G R (z=0, k). It is easy to find if one looks at Eq. (12. 6 j) at t=0, 

POD 

(n(t=0, k)> = -/io(k) / dt' e^'G R n (t', k) = -^ (k)G« n (^=0, k) . 

Jo 

So in the small-k limit we can identify G R n (z=0, k) = — %, which gives the retarded function 



CK^k) 



This function is analytic in the upper-half plane of complex z as it should be. We can 
define G r (lu, k) in the whole complex plane as the analytic continuation of G R (z, k) from the 
upper-half plane. In the lower half-plane, the retarded function has a pole at u — —iDk 2 , 
corresponding to the diffusive mode. 

Given the above expression for the retarded function, we can deduce the standard Kubo 
formula for the diffusion constant, 

Dx = - lim lim — Im G R n {u, k) . 

lj^O fc-s>0 k 

2.2 Canonical approach to hydrodynamic response functions 

The above example of diffusion allows for a simple generalization when there are several fields 
present. Let <p a (t,x) be the set of hydrodynamic variables which have microscopic operator 
definitions, such as the charge density J°, momentum density T° l etc. We add the sources A a 
to the Hamiltonian as 

SH = - I d d x A a (t,x)^ a (t,x). 
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The master relation (12 .4p now gives 

<%> (t,x)> = - Hdi' jd d x' GK(t-t',X-x')\ b (t',x'), 

J — oo ./ 

or in terms of the Fourier components simply 

6(<p a (u, k)> = -G«(u, k) X b (u, k) . (2.8) 

As before, the retarded function is defined by 

G*{t-t, x-x') = -i9(t-t') (& a (t, x), <p b (t', x')]> . (2.9) 

We can repeat the initial- value analysis just like we did for diffusion. In momentum space 
the variables ip a obey linear equations 

dt(p a (t, k) + M ab {k) Vb {t, k) = , (2.10) 

where the matrix M ab is determined by the equations that ip a satisfy, such as the equations 
of relativistic hydrodynamics. Applying the Laplace transform, one finds 

(-izSab + M ab )(p b (z, k) = (p° a (k) , 

where y2°(k) = tp a (t=0, k). The initial values of the hydrodynamic variables are related to 
the initial values of the sources by the susceptibility matrix, y3°(k— >0) = \ab ^°(k— >0), wher^f] 



Xab 



d(p a 



is the static thermodynamic susceptibility. As before, the retarded function at zero frequency 
and small momentum is just minus the susceptibility, G^ b (z—0, k— >0) = —Xab- The solution 
to the hydrodynamic equations can now be written as 

Va (z,k) = (K- 1 ) abXbc X° c (k), (2.11) 

where K ab = —iz5 ab + M ab (k). On the other hand, <p a ( z , k) can be expressed in terms of the 
retarded function, just like we did for diffusion: 

<p a (z, k) = -- (G&(z, k) - GS(z=0, k)) Ag(k) . (2.12) 

iz v ' 

Comparing the solution to the hydrodynamic equations (12.111) with the expression (I2.12p . we 
find the retarded function 

G R (z,k) = -(l + tzK- 1 )xi (2.13) 

omitting the matrix indices. 



5 Note that the relation between the source and the field is local in space, but non-local in time, so it would 
be incorrect to write </ 3 a(w, k) = Xab Af,(a;, k). 
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2.3 General properties of response functions 

In addition to the retarded function (I2.9p . other response functions are useful. For bosonic 
operators ip a , <Pb, we define 

G£(t-f ,x-xO = -i6(t-t') ([^ a (t,x),^(t',x')]> , 

G^(t-t',X-x') EE i0(f-t) ([^ a (t,x),^(t',x')]> , 
G a6 (t-t', X-X') EE i ({^ a (t, X), <p b (t', x')}> , 
p afe (t-t',X-x) EE ([^ a (t,x),^ b (t',x')]) . 

The first one is the retarded function, the second one is the advanced function, the third one 
is the symmetrized function. The Fourier transform of the commutator p a b{uj, k) is called 
the spectral function. The expectation values are taken in static thermal equilibrium in the 
grand canonical ensemble, and the Heisenberg operators are defined with the Hamiltonian 
H' = H — fiQ. The above four functions are not independent, and can all be deduced from 
the Euclidean time-ordered function, see e.g. [38J. 

Thanks to the theta functions of time, G R (u, k) is analytic in the upper half plane of 
complex co, while G a (oj, k) is analytic in the lower half plane of complex uj. Translation 
invariance implies that 

G«( W ,k) = Gi(-w,-k), 
pab(u, k) = -p ba (-u, -k) . 
From the definition of G R and G one finds 



n R,A( in f du} 



2tt uj — uj' ± is ' 

with the upper sign for G R and the lower sign for G . For Hermitian operators (p a and (pb, 
the matrix p a b(u),k) is Hermitian. The diagonal components p aa (uj,\i) are then real, and we 
find for real uj 

ReG^u, k) = ReGf>, k) = P f P -^l , 

J Z7T UJ — UJ 

lm G R a (uj,k) = -ImG&(w,k) = -\p aa (uj,\s) , 

where P denotes the principal value of the integral. 

Writing out the definitions of G a b(u>, k) and p a b(u, k) in the basis of H' eigenstates in the 
grand canonical ensemble and inserting a complete set of states, one can easily show that 

1 i+ e _/3w 
G ab (uj, k) = - y-— ^ p ab (u, k) . 

In the hydrodynamic regime (3uj <C 1 this gives p a b(uJ, k) = (3ujG a b(uj, k), or 

IT 

G ab (uj,k) = ImG*(u;,k). (2.14) 

UJ 



20 



For a Hermitian operator (p a (t,x), the spectral decomposition implies G aa (u;,k) ^ 0, or 
^Paa(w,k) ^ 0, which gives 

-ImG^(w,k) ^ for u ^ 0. (2.15) 

Applied to hydrodynamic response functions, this condition will imply that the transport 
coefficients have a definite sign: 77 ^ 0, £ ^ 0, and cr ^ 0. 

In addition, the retarded functions have to be consistent with the symmetries of the 
theory, such rotation invariance, parity, charge conjugation, time-reversal, and various global 
symmetries. Time-reversal turns out to be particularly useful in constraining the transport 
coefficients that appear in the constitutive relations, as was pointed out by Onsager [39J SO] . 
Hence the constraints on the transport coefficients due to the time-reversal covariance are 
called "Onsager relations" . The implications of time-reversal for the retarded function ( 12.91) 
are easy to derive. The ant i- unit arity of the time-reversal operator O implies that jH] 

(P\A\a) = (aie^e- 1 !^) 

for any linear operator A, where \a) = Q\a) are the time-reversed states. We now apply 
this to the retarded function (12.91) of Hermitian operators ip a which transform in a definite 
way under time-reversal, 0<^ a (t,x)0 _1 = rj a <p a (—t,x), where r] a = ±1 is the time-reversal 
eigenvalue of ip a . The expectation value in (12. 9ft is taken in the grand canonical ensemble, 

(...) = -tr(e-^ +/3 ^...), 

and the trace can be taken either in the basis of energy eigenstates |n), or in the basis of the 
time-reversed states \n). If the microscopic system is time-reversal invariant, i.e. [H, 0] = 0, 
then the space-time translation invariance implies that 

G%{t, X) = -ld{t) {[(p a (t, X), <p b (0)]) = -i0(t) (Mt, -X), <p a (0)]) T) a r) b = GH(t, -X^aTfc . 

If the microscopic system does not have time-reversal invariance, it may still be possible to 
use time-reversal covariance to find a similar relation for the retarded function. Namely, let 
us assume that the Hamiltonian depends on some time-reversal breaking parameters B such 
that QH(B)Q~ 1 = H(—B). For example, B could be the external magnetic field, or it could 
be the mass for Dirac fermions in 2+1 dimensions. The retarded function must then satisfy 

G£(w, k; B) = ThVtG^u, -k; -B) . (2.16) 

This can be also written in the matrix form as G(u, k; B) = SG T (u, — k; —B)S , where S = 
diag(?7i, 772, • • • ) is the matrix of time-reversal eigenvalues of the hydro variables (f a , which 
satisfies S 2 =l. Equation (12. 16[) is the basis for the Onsager relations. Let us now apply it to 
our hydrodynamic correlation functions. 

Taking u=0 and k— s-0 in (" 12 . 1 6[) we find for the static susceptibility matrix 

S X (B) S = X T (-B) . (2.17) 

Further, applying the relation (I2.16P to the hydrodynamic correlation function (I2.13p . evalu- 
ated in the linear response theory, one finds 

X (B)SM T (-k; -B) = M(k; B) X (B)S . (2.18) 
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In these lectures, we will not consider any T-breaking parameters B; in fact the breaking 
of T due to either magnetic field, or the fermion mass in 2+1 dimensions would give rise to 
extra terms in the constitutive relations (|1.16p . see Ref. [16J. The condition ( 12.18)) may be 
viewed as a constraint on hydrodynamic constitutive relations arising from the time-reversal 
covariance. Since the matrix M(k) simply represents the constitutive relations in linearized 
hydrodynamics, this tells us that the constitutive relations in hydrodynamics can not be 
arbitrary, but must be such that Eq. (" 12 . 1 8[) is satisfied. As an example, we will see shortly 
that the constitutive relations (I1.16P in relativistic hydrodynamics do not satisfy Eq. ( I2.18p . 
unless Xt — 0. In other words, time-reversal covariance may be used to deduce that Xt — 0, 
without alluding to the entropy current argument. 



2.4 Retarded functions in relativistic hydrodynamics at [i = 

As a simple example of the above general formalism, let us consider small hydrodynamic 
fluctuations about an equilibrium state with ft = 0, in other words n = in equilibrium^ 
As the Eckart frame is ill-defined for states with n = 0, we will use the Landau frame. The 
hydrodynamic equations are thus given by the conservation laws ( II. 7p . supplemented by the 
constitutive relations (I1.16p . We linearize the equations around the static equilibrium state 
v l = 0, T = const, /i = 0, and choose the hydrodynamic variables to be the fluctuation in 
the energy density 5e(t,x) = 5T 00 , momentum density iii(t, x) = T°\ and charge density 
n(t,x) = J°. 

First, let us look at the constitutive relation for the current in (I1.16p . All terms in the 
constitutive relation must transform the same way under charge conjugation C as the current 
J M does, in other words a is C-even, and Xt is C-odd. Assuming a C-invariant microscopic 
theory, the only source of C violation in hydrodynamics is the chemical potential /i, hence 
Xt — > as /i — > 0. In other words, when /i vanishes in equilibrium, the Xt term does 
not contribute to the linearized hydrodynamic equations. The fluctuation in charge density 
n(t, x) then decouples from the fluctuations of energy and momentum densities, and satisfies 
the diffusion equation which we already studied in Section 12.11 The diffusion constant is 
D = cr/x, where x — (dn/dfj,)^ is the static charge susceptibility. Hence, the two-point 
retarded function of J°(t, x) is given by 

The equations for Se(t, x) and 7Tj(i, x) come from the conservation of the energy-momentum 
tensor. We Fourier transform in space, and decompose 7Tj = 7rJ + irj-, where parallel and 
perpendicular are with respect to the direction of the spatial momentum. Choosing k along 



6 As a reminder, a state in the canonical ensemble which has a fixed value of the total conserved charge 
(such as the number of particles) is not the same as the state in the grand canonical ensemble at /i=0. 
Hydrodynamic fluctuations in a state with n — is a perfectly sensible thing in relativistic hydrodynamics: 
n = simply means that there is an equal number of particles and antiparticles in thermal equilibrium. On 
the other hand, in non-relativistic systems there are no antiparticles, and n = means that the equilibrium 
state is empty, with no matter that could flow. 
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x, we have the following set of equations 

d t Se + ik x iT x = , 

<9t7r" + ik x v 2 s 5e + 'j s k 2 7f x = , 

d&i + 7r,kV = • 

Here v 2 = dp/de, 7^ = T]/(e+p), 7^ = (^%-^^+C)/(^+p); V an d C are shear and bulk viscosities, 
and d is the number of spatial dimensions. 

The transverse momentum density 717- obeys the diffusion equation, with 7^ playing the 
role of the diffusion constant. This means that we should be able to use the results for the 
diffusion equation. The source for the momentum density is the velocity, and it appears in 
the Hamiltonian in the form 



5H = — J d d x Vi(t, x) 7Tj(t, x) 



[think grand canonical density operator, —/3H — > fiu^P^ ~ —/3(H—\-P)]. The role of the 
susceptibility x is played by the enthalpy density w = (e+p), because to linear order 7Tj = wvi. 
The retarded correlation function of the transverse momentum density is now easy to write 
down: for k along x, we have 

G^>,k) = 7 ^^. (2.20) 



icu — 7 r) k 2 



Now let us look at the coupled equations for 5e and 7P . These equations describe sound waves, 
and have the form of Eq. (12.101) with ip a = (5e, ir x ), and 

( ik x 

M a b = I 2 2 

V«Abv; 7* k 

To find the susceptibility matrix, we need to identify the sources A a . A disturbance in energy 
density can be created by a disturbance in temperature, —f3H — >• — (/3 + 5(3)H. Noting that 
6p/p = -8T/T, we have 

6H = - jd d x ^ r ^' X ^ e(t,x) +u i (t,x)7r i ft,x)^ . 

Thus the sources corresponding to the fields y? a = (5e, 7r x ) are A a = (8T/T,v x ). The suscepti- 
bility matrix is therefore diagonal, 



Xa6 



c^T 
u) 



where c v = de/dT = (H 2 ) conn /(VT 2 ) is the specific heat. At /i = 0, it can be easily related 
to the enthalpy and the speed of sound, w/T = s = ^ = ^-§f = v 2 c v . The matrix of 
T-eigenvalues is S — diag(l, —1), and one can easily check that the relation (" 12 . 1 8[) is satisfied, 
so the response functions will come out consistent with the time-reversal invariance. The 
retarded function is given by Eq. (I2.13p . 

r,R ( i\ ™ f ^ uk * 

G ab {u, k) - 



to 2 — k 2 t> 2 + iuj'jsk 2 \u)k x k 2 v 2 — ioujsk 2 
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In the plane of complex oj, there are two poles, when the denominator vanishes. In the limit 
of small momenta, the poles are at u> = ±|k|u s — i'j s k 2 /2, corresponding to weakly damped 
sound waves. We can combine the contributions from n 1 - and 7r" into one retarded function 
for momentum density, 

g?.,.(u,, k) = ( s iS - k M + f k !5 :T 7 ' k l ■ < 2 - 2i > 

^ v ; \ J k 2 J iuj - 7,,k 2 k 2 w 2 - k 2 w 2 + z^7 s k 2 
For completeness, the other retarded functions are 

Gf^,k) = G*>>k) * uki 



uS 1 — k 2 f 2 + zco>7 s k 2 
w k 2 



w 2 — k 2 f 2 + ia;7 s k 2 

Evaluating the imaginary parts of the retarded functions, we find the Kubo formulas 

a = -^ImGi(u;,k^0) , (2.22a) 

1 = -^p (*« - ^r) k^O) , (2.22b) 

—n + C = - F Im G£(w, k^O) . (2.22c) 

The positivity condition (12.151) implies that 77 ^ 0, £ ^ 0, and a ^ 0, consistent with the 
requirement that small hydrodynamic fluctuations decay (rather than grow) with time. Using 
the relation ( 12. 14(1 . the above Kubo formulas can be equivalently expressed in terms of either 
the symmetrized correlation functions G a b, or the spectral functions p a b- 

The Kubo formulas can be also written in terms of correlation functions of spatial currents, 
rather than charge densities. To do so, we demand that the correlation functions of conserved 
currents satisfy 

PG JflJ >,k)=0, 
FG^( u ,k)=0. 

For the conductivity, we then find 

° = Wd Gj ^ {u ' k=0) = 2Y Gj ^ (a; ' k=0) • (2 ' 23a) 
In writing down this expression, we have used the relation ( 12. 14(1 between the symmetrized 
and retarded functions, as well as time-reversal invariance which by Eq. (12. 16(1 implies Gj j i = 
Gjjjo, and rotation invariance which implies that at zero spatial momentum Gj i j j must be 
proportional to Sij. A similar argument gives Kubo formulas for shear viscosity 

77 = — H™{ G T . T . = — G T T , (2.23b) 



and for the bulk viscosity 
$ S G 



(2.23c) 
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where = + \5\5l, — \^ l ^ki is a projector onto symmetric traceless tensors, and 

all correlation functions are evaluated at k=0. The positivity of bulk viscosity implies that 
zero-momentum correlation functions of the stress tensor obey 

2d-2 

G TxxTxx (u,k=0) > —G TxyTxy (u,k=0), 

as a consequence of rotation invariance. Taking a closer look at the hydrodynamic retarded 
functions, we see that 

uG%u,k)-k j G« e (u,k) = 0, 

wG* .(w, k) - kjG^.(u 7 k) = wh , 

in other words there is a contact term in the right-hand side of momentum conservation 
equation applied to our correlation functions. 
Finally, let us look at the the speed of sound: 

s {ds/dT) ' 

as follows from the relation s = v 2 s c v . While it is easy to see that v 2 ^ because both s 
and c v are non- negative, thermodynamic inequalities do not demand that v 2 ^ 1, as one may 
expect in a consistent relativistic theory. Demanding that v 1 ^ 1 produces a constraint on 
the equation of state. For example, for the equation of state of the form p(T) = cT a , one 
must have a > 2 in order for the speed of sound to be less than the speed of light 

2.5 Retarded functions in relativistic hydrodynamics at \x ^ 

Now let us consider linearized relativistic hydrodynamics at non-zero \x in equilibrium. This 
means that the equilibrium state has non- zero charge density n, and as a result, the fluc- 
tuations of charge density will couple to the fluctuations of energy and momentum density. 
Again, the hydrodynamic equations are given by the conservation laws (11. 7J) , supplemented by 
the constitutive relations (|1.16p . in the Landau frame. We linearize the equations around the 
static equilibrium state v 1 = 0, T = const, /i = const, and choose the hydrodynamic variables 
to be the fluctuation in the energy density 5e(t, x) = 5T 00 , momentum density 7Tj(t, x) = T 0t , 
and charge density n(t, x) = J°. Choosing the spatial momentum in the x direction, the 
conservation of the energy-momentum tensor gives 

d t 5e + ik x 7i x = , 

d t 7i^ + ik x Pi5e + ik x (3 2 5n + j s k 2 ir x = , 

where Pi = (fjr) , P2 = (f^) ■ The equation of current conservation gives 

dt5n + —ik x n x + o~a:ik 2 <5e + aa2k 2 5n = , 

w 

7 More precisely, the effective "speed of light" which determines the Lorentz invariance of the microscopic 
theory. For example, the effective "speed of light" in graphene is about 300 times less than the velocity of 
electromagnetic waves in the vacuum |42j . 
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where the coefficients a\ and ot2 are 



djj, 
d7 



V Xt 
T a 



«2 



<9n 



^ Xt 

T o 



dn J 



The equation for n decouples, as before, and gives rise to the same two-point function ( I2.20p . 
The remaining equations have the form of Eq. ( I2.1(jp with (p a = (Se, u x , Sn), and 



M, 



ik x 
„/, = | ik x (3x 7 s k 2 //.•,. I, 
^octik 2 ik x n/w aa^k 2 



Again, to find the susceptibility matrix, we need to identify the sources A Q . For infinitesimal 
constant disturbances ST, Sfi, and Vi, the Hamiltonian changes as 



H -^H 



5T 
~T 



H-fiQ) -SfiQ - v-P 



For non-constant slowly varying sources, we take 
fST(t,x 



SH 



/a ( £(t , x) - „ B( «, x) ) + Mt , x) n{t , x) + „ (t , « X ) 



Thus we identify the sources corresponding to (p a = (Se, ir x , Sn) as A a = (ST /T, v x , Sfj, — £ST), 
and therefore the susceptibility matrix is 



Xab 



w 



dfi ) T 




rp I dn \ 
1 VdT/ 



H/T 







dn 

8fJ, J rp 



(2.24) 



Note that T(^)^ T = T(^)^ + /^fj^l > and similarly for n. The definition of the thermo- 
dynamic quantities in the grand canonical ensemble implies that T(dn/ cfT)^/? = (de/dfi)x, 
in other words, the susceptibility matrix is symmetric, X13 — X3i> consistent with (I2.17p . 
Further, Xu ^ and X33 ^ because the first is proportional to (H 2 ) conn , while the 
latter is proportional to (iV 2 ) conn . In addition, det(x) ^ because it is proportional to 
(H 2 ) conn (N 2 ) conn — (HN)l onn , which is non-negative by the Schwarz inequality. Thus the 
susceptibility matrix \ab is positive-definitejf] 

An exercise in thermodynamic derivatives shows that 

#1X11 + A2X31 = w, 
AX13 + P2X33 = n, 
ociXu + «2X3i = -Txt/ v , 

«lXl3 + "2X33 = 1 • 



8 Note that only those p(T, /i) that follow from the partition function in the grand canonical ensemble 
represent legitimate equations of state. An arbitrary function p(T, fi) will give rise to a susceptibility matrix 
which is not necessarily positive-definite. 
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The matrix of T-eigenvalues is S = diag(l, —1,1), and it's a matter of simple algebra to check 
that the condition of time- reversal covariance (12 . 1 8[) can only be satisfied if Xt = 0, which 
we take from now on. The coefficients 0i i2 and ai$ can be easily expressed in terms of the 
components of the susceptibility matrix Xab- 

The retarded functions can be found by Eq. (12. 13ft ; in the limit u— >0 and k— >0 one finds 

vo k 2 

G ^' k) = ^k)^ + ^ 2k2 )' 

k 2 / \ 

G nni^^) = d ^ k ^ (un 2 /w + iak 2 (a 2 l3 1 -a 1 l3 2 )x33j , 

GiM k) = k) = ^ + ^aa^) , 

G* xn (u, k) = G^(u, k) = — ^ (nu - ifxaiwk 2 ) , 

k 2 

G? n (u, k) = G£(w, k) = (nuj - za^wk 2 ) . 

where 

d(u, k) = u 3 + zw 2 k 2 (7,+aa 2 ) - wk 2 w 2 + 2ak 4 (ai/3 2 -a 2 /3i) , (2.25) 

and v 2 = /?i + n(3 2 /w. Setting (i(u;, k) = determines the eigenfrequencies of the system. 
There are three modes, whose frequencies in the limit k — > are 

u = ±kv s ~ i — k 2 , 



-iDk 2 



where T = 7 S + a(3 2 («i + ^2) /v 2 , D = a (a 2 f3i—ai(3 2 )/v 2 . These are the familiar sound 
and diffusive modes. The expression for the speed of sound in terms of the components of the 
susceptibility matrix (I2.24p is 



2 = n 2 xn + w 2 X33 ~ 2rawxi3 
det(x) 

The numerator can be written as (n^/xn ~~ w^/xzz) 2 + 2nu)( v /xnX33 — X13)) an d therefore 
the positive-definiteness of the susceptibility matrix implies v 2 ^ 0. By the same argument 
as in Sec. 11.21 one finds v 2 s = (dP/de) s N . Note that the speed of sound does not depend on 
the dissipative transport coefficients, and is only determined by the equation of state. In a 
scale-invariant theory, the equation of state p(T,fi) = T d+1 f(T/n) implies p = e/d, which is 
equivalent to the tracelessness of the energy-momentum tensor (11. lip in equilibrium. Thus in 
a scale- invariant theory, the speed of sound is v s = l/y/d, and is independent of temperature 
and chemical potentials The expression for the diffusion constant in terms of the components 



9 In QCD, the equation of state at asymptotically high temperatures is p — e/3, and therefore the speed of 
sound at asymptotically high temperatures is 1/ \/3- 
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of the susceptibility matrix is 



D 



aw 2 



n 2 xn + w 2 X33 ~ 2nwxi3 ' 



and shows that D is only positive for a positive. The expression for the sound wave damping 
constant in terms of the components of the susceptibility matrix is 

r= 1 ( 2d ~ 2 r] i 1 a ™ (^xii-w%3) 2 



w \ d J det(x) n 2 xn + w 2 X33 - 2rai>Xi 3 

Evaluating the imaginary parts of the retarded functions, we find the Kubo formulas (I2.22p 
for the transport coefficients. In other words, the Kubo formulas for rj, (, and a have the 
same form at fi 7^ as they do at fx = 0, even though the values of the transport coefficients 
of course depend on fi. Again, the positivity condition (I2.15P implies that 77 ^ 0, ( ^ 0, and 
a ^ 0, consistent with the requirement that small hydrodynamic fluctuations decay (rather 
than grow) with time. Our hydrodynamic correlation functions satisfy 



< - = , (2.26a) 

wG>, k) - kjG^u, k) = , (2.26b) 
k) - kjG^u, k) = wk t , (2.26c) 



again with a contact term for momentum conservation in the right-hand side. 

Similar to the sound in the fi = fluid, thermodynamics does not constrain the speed of 
sound to be less than the speed of light. Again, demanding that v 2 ^ 1 produces a constraint 
on the equation of state p(T, fx), 

n 2 Xn + w 2 X33 ~ 2mZ?xi3 < det(x) • 



2.6 Variational approach to hydrodynamic response functions 

The response functions we computed so far were evaluated using the standard linear response 
approach, as described for example in Ref. [2S]. After introducing external sources for the 
conserved densities (energy density, momentum density, charge density), the response func- 
tions can be evaluated starting from Eq. (12 .4p in the canonical (operator) formalism. The 
advantage of this approach is that the external sources can be introduced quite easily: their 
coupling to conserved densities follows from equilibrium thermodynamics. The disadvantage 
is that the response functions of spatial currents (charge current, stress tensor) are not eas- 
ily accessible. For example, the component of the spatial current which is transverse to the 
spatial momentum does not couple to the charge density, and hence the correlation function 
(8ij — kikj/k 2 )G jijjioj, k) at non-zero k can not be determined by this method. 

The general relation (12. 8p between the fields and the sources suggests that it would be 
more elegant to introduce sources for J M and rather than just for the conserved densities. 
The retarded functions G^ h {u, k) can then be evaluated by taking the variation of one-point 
functions with respect to the source. The advantage of this approach is that it gives direct 
access to all response functions G j j , Gj T , and G? T . The disadvantage is that it may 
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not always be obvious how to couple the currents to external sources. In relativistic hydro- 
dynamics, Lorentz symmetry is of great help. In what follows, we will write the constitutive 
relations in the presence of background metric and background gauge field which cou- 
ples to the conserved current J M . When u^, T, and \i satisfy the hydrodynamic equations, 
the constitutive relations provide the one-point functions of and in the presence of 
external sources. Defining 



JHx) 



-g (J^x)) 



V 



pi/ , 



X 



-g (T^(x)) Ag . 



(2.27) 



and expanding the metric about the flat space as g^ v = f]^ u +h^ u , the retarded functions may 
be defined as 



G jarp^v ( X J 



-2 



6A„(p) 
bjHx 



■fW, 



X) 



A=h=0 



<*V(o) 



^R 

Tot 



r p<3"T r pi- Ll/ 



X 



-2 



SA a (0) 
5T aT (x) 



A=h=0 



A=h=0 



(2.28a) 
(2.28b) 



A=h=0 



The factor of \f—g gives rise to contact terms, by virtue of S^—g 



^gg^Sg^. In the 



Landau frame, the constitutive relations in the presence of external sources are 

rpflU 



eu>V + pA^ - 77 A> J,a A v P (v a up + V p u a - - d gapV^J - CV M w" , 



J" 



nw 



+ aA^Vx + X^ X Ex + XtA^VxT 



(2.29a) 
(2.29b) 



where A^ = g^+u^u", V a is the covariant derivative, V x = E x -TA Xp V p (n/T), = F pu u v 
is the electric field, and is the field strength of the external gauge field. None of the 
hydrodynamic variables carry charge, hence there are no gauge-covariant derivatives. The 
velocity satisfies = — 1, and consequently u a V M u" = because the metric is covariantly 
constant. In writing down the constitutive relations for the current we have grouped the 
three independent vectors VaT, Va/U and E\ into combinations with coefficients a, Xt, and 
Xe- The coefficient a is the electrical conductivity, and the term aA^Ex is just a relativistic 
version of the Ohm's law. As for the coefficients Xt and Xe (which are allowed by Lorentz 
covariance), they must vanish as a consequence of either the positivity of entropy production, 
or of consistent thermodynamics with external sources, see for example [16, [36] . 

With the external sources A p and g^, the hydrodynamic equations take the form 



F vX J x 







(2.30a) 
(2.30b) 



In order to find the retarded functions, we linearize the hydrodynamic equations around the 
equilibrium state with constant e, p, and n, and solve the conservation equations (12.30P with 
the constitutive relations (|2.29|) . expressing ST, dp,, and Su^ in terms of A^ and h pu to linear 
order. Upon substituting the solution back into (I2.29p . the retarded functions can be read off 
from the definition (12. 28 p . 

The retarded functions obtained in this variational approach will differ from their coun- 
terparts in Section 12.51 by contact terms. For example, taking the spatial momentum in the 
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x-direction, we have for the transverse momentum and current densities 

7/k 2 
iuj — 7^k 2 



G T o vT oy (u>, k) = — + e , (2-31) 



n 2 iuj 



G% Jy (u, k) = -iua + — — -. — , (2.32) 

e + p iuj — 7^k 2 

where 7^ = 77/(5 + as before. Note that at non-zero charge density in equilibrium, Gj yJy 
has the same shear-mode singularity as G^, 0yTOy , due to the convective term in the current. 
For the retarded functions of spatial currents at zero spatial momentum we find 

ft 2 

G% JX = iua + . . . , (2.33a) 

e + p v ; 

G^x yT x y = p — icorj + . . . , (2.33b) 

1 J) 

- (Gtxx T xx + (cI—1)Gt*xtvv) — ~^ + + P) v s ~ + • • • • (2.33c) 

Remembering the relation f)2.14p between G a b and G^ b , this gives precisely the Kubo formulas 
(I2.23p . The variational response functions satisfy 

cuGj^oo jo — kjGrpoj jo = , (2.34a) 

ujGj*qoj<qo — kjGj^ojj^oo — — ecu, (2.34b) 

oj&poorpoi — kjGj-ojrpQi = p fcj , (2.34c) 

with contact terms in the right-hand side. Comparing with the analogous expressions (12. 26ft . 
we see that the contact terms differ between the response functions evaluated in the canonical 
approach and in the variational approach. For a discussion of contact terms in the variational 
approach, see Ref. [43J. 

We conclude the section about correlation functions with a comment about the choice 
of variables in the constitutive relations. When writing down the constitutive relations in 
( 11.16P or ( 12.29(1 . we chose to use the transverse vectors A mA VaT and A mA Va/U. However, one 
could have chosen to work with a different set of transverse vectors, say with A mA VaT and 
A^iix, which is a choice adopted by Eckart [33], [28] - By using the zeroth-order hydrodynamic 
equations, we find 



E» - TA^ a Va(|)] = ° € Af^ X [T^x + V A T] + 0{d 2 



The latter form of j M is, however, not suitable for studying hydrodynamic fluctuations in the 
"neutral" state with n = 0, hence our preference for working with A^ a VaT and A^Vxfi. 

Another place where this choice of variables in the constitutive relations becomes peculiar 
are response functions in a state with non-zero n, evaluated from Eckart-frame hydrody- 
namics. Focussing on the response function of the transverse momentum density (the shear 
mode), and using Eckart-frame constitutive relations (11.1 Tj) with g M = —cxV fJ- (e+p)/n leads 
to expression ( 12 .3 1(1 . as before. On the other hand, using Eckart-frame constitutive relations 
(Tl~iTj) with q» = -kA» x [Tii x + V A T] leads to 

G^ 0yT o y (uj, k) = . j- -. — . ' J °^ i 2 + ^ ' (2.35) 
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where uq = n 2 /(aw), and we have used the relation between heat and charge conductivities 
K = a(e+p) 2 /(n 2 T). There is now a 0(u 2 ) term in the denominator of the response function, 
however, in the hydrodynamic regime, u/coo <C 1, and this term should be ignored^ If one 
naively proceeds with the analysis of the singularities of the response function (12.35p . one 
finds two poles: the normal shear pole at lo\ = —i^k 2 + 0(k 4 ), and the second pole at 
U2 = icoo + 0(k 2 ). The second pole is in the upper half-plane of complex u>, and, if taken 
seriously, would represent an instability of the thermal equilibrium state [37J. The instability 
is, however, fictitious, as the second pole is outside of the validity regime of hydrodynamics: 
the corresponding time scale a(e+p)/n 2 is microscopic, not macroscopic (in a gas, this time 
scale would be determined by the mean-free time between particle collisions). Hydrodynamics 
is only a sensible effective description at time scales much longer than such microscopic times, 
and the correct way to interpret the response function (12.35p is to expand (u— 0J2) in the 
denominator in powers of uj. This gives 

Gfo aT o>, k) = ^ ff 2 +Q(k4) I 1 + ( k2 ) + °^ 2 ) + °^ 2 )} + ^ 

and shows that the contribution of the term u/uo in ( 12.351) is of the same order as contri- 
butions from higher-order hydrodynamics. Therefore, this term should not be retained in 
response functions evaluated using first-order hydrodynamics. In first-order hydrodynamics, 
the response function of the transverse momentum density is given by Eq. (I2.3ip . and only 
has a pole in the lower half-plane of complex u, as it should. It happens quite generally 
that certain terms in response functions evaluated from first-order hydrodynamics may be 
"contaminated" by higher-derivative corrections. Thus, some care needs to be taken in order 
to determine the response functions reliably, up to a given order in the derivative expansion. 



3 Interactions of hydrodynamic modes 

In the last section, we looked at two-point correlation functions of T^ v and J M in thermal equi- 
librium state, evaluated from linearized hydrodynamics. What happens to these correlation 
functions if the non-linearities in the hydrodynamic equations are taken into account? While 
the linear equations describe freely propagating hydrodynamic modes, the non-linearities are 
responsible for the interactions among the modes. These interactions will modify the equi- 
librium correlation functions of T^ v and J M , and in particular will modify the transport 
coefficients 77, (, and a compared to their naive "bare" values in linearized hydrodynamics. 

The modification of viscosity compared to its naive value can be understood as follows. 
An inhomogeneous velocity profile v x (y) of the shear mode will eventually equilibrate due to 
particles transferring rc-momentum between the layers of fluid, as shown in Fig. [TJ However, 
momentum can also be transferred by collective excitations such as sound and shear modes, 
with wavelength shorter than the length scale of the inhomogeneity of the flow, as illustrated 
in Fig. |2j Such collective excitations will be present even in thermal equilibrium due to fluctu- 
ations, hence the shear viscosity measured for example from equilibrium two-point correlation 
functions will receive contributions from the collective modes. To estimate the contribution to 

10 Second-order terms in the constitutive relations will only contribute 0(d 3 ) terms such as 0(wk 2 ), 0(u> 3 ) 
to the denominator of (|2.35[l . 
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Figure 2: A fluid flow with an inhomogeneous velocity profile v x (y) corresponding to the 
shear mode, represented by thick arrows. Thin wiggly arrows represent sound and shear 
modes generated by thermal fluctuations in fluid elements in local thermal equilibrium. These 
collective excitations can transfer x-momentum in the ^-direction, thus contributing to shear 
viscosity. 

r\ from the sound waves, one can think about the short-wavelength sound waves as particles, 
whose "mean free path" (propagation distance) is proportional to l/(-y s k 2 ), as is evident from 
the sound wave dispersion relation u> = ±kv s — i , ~f s k 2 /2. Hence the contribution to 77 from 
sound waves is proportional to J d d k/ '( , y s k 2 ), in other words is inversely proportional to rj it- 
self. The momentum integral is UV divergent at large k and needs to be cut off at some fc max , 
which is the maximum value of k at which the hydrodynamic picture of collective excitations 
is still applicable. 

In this section, we will give a simple derivation of how two-point correlation functions 
of conserved currents are modified by the non-linear terms in the hydrodynamic equations, 
and estimate the corresponding fluctuation correction to the shear viscosity. We will also 
see that hydrodynamic fluctuations will lead to the breakdown of the derivative expansion 
in hydrodynamics. Put differently, second-order hydrodynamics (taken as a set of classical 
partial differential equations) does not correctly describe the relaxation of T^ v and J^ 1 in the 
hydrodynamic regime. In the context of non-relativistic hydrodynamics, the problem with 
the derivative expansion was known for a long time [23], and relativistic hydrodynamics is no 
different in this respect. It is important to emphasize that the physics of hydrodynamic fluc- 
tuations is only partially contained in the equations of hydrodynamics if the latter are taken 
merely as classical partial differential equations. This is because the classical equations of hy- 
drodynamics ignore the excitation of collective modes by thermal fluctuations: the classical 
equations contain only the "dissipation" part of the fluctuation-dissipation theorem (Hj. 

3.1 Correlation functions of conserved currents 

Following Ref. [I5|l2l]. we will evaluate the effect of the non-linear terms in the hydrodynamic 
equations on the correlation functions of spatial currents J 1 and stresses T 4J in the hydrody- 
namic regime u — > and k — > 0. For simplicity, we will do this in a scale- invariant theory, 
in the equilibrium state with [i = 0. Choosing the hydrodynamic variables as 6e = T 00 — e, 
7T* = T° l , and n = J°, the constitutive relations in the Landau frame take the following form 
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to quadratic order in fluctuations: 

T ij =(p + v 2 s 5e) 5 ij - 7„(<9V + dV - § 6 ij d k ir k ) + 4#2 vrV + ... , 

f = -D &U + — + . . . . 

w 

Here w = (e + p) is the equilibrium enthalpy density, D = a/x is the diffusion constant for 
charge density, and j v — 7)/w is the diffusion constant for the transverse momentum density. 
The speed of sound is v 2 = 1/d, and the bulk viscosity vanishes due to scale invariance. The 
projector onto traceless symmetric tensors is H 1 ^ = + h^i^i ~ h^Ski- The dots denote 
the terms which are suppressed either due to containing higher powers of the hydrodynamic 
variables, or due to containing extra derivatives of the hydrodynamic variables. The contri- 
butions to the equilibrium correlation functions from the quadratic terms in the constitutive 
relations thus have the form 

G% k (t, x) = ij(n(f , x) tt% x) n(0) n k (0)) , 
G% Tki (t,x) = -^^ n ^(7r-(t,x)7r-(^x)^(0)^(0)). 

Assuming that the small equilibrium fluctuations are Gaussian and therefore can be factorized, 
we have for the quadratic contributions to the connected correlation functions 

GS>.k) = ±J ^I^GU^kOG^-a/.k-k'), 

= h H>L < I ^Sr G °'-- (uj '' k ' )G: -*>- w ' k " k,) ■ 

Here G^ n (u, k) and G° (u, k) are the symmetrized (or unordered) correlation functions, 
obtained in the linear response theory. Explicitly, relations f )2.19p . f)2.2ip . and (I2.14p give 



2TD X k 2 

™ ( ' } = u 2 + (Dk 2 ^ 2 ' ( } 



' J v k 2 y u 2 + ( 7 ,k 2 ) 2 + k 2 (u 2 - k 2 v 2 ) 2 + (w 7s k 2 ) 2 ■ 1 } 

Let us first look at the current-current correlation function Gj)j (u, k). There are two contri- 
butions, one of the "shear-shear" type, and one of the "shear-sound" type. After performing 
the frequency integrals, one is left with momentum integrals which are linearly divergent in 
d — 3, and logarithmically divergent in d = 2. Introducing the high-momentum cutoff A, one 
finds at zero momentum and small frequency in d = 3 



v m nm a 2 (2T x )(2T^) H 1 / 2 

G ^>' k=0) = ° (A) " % 4^ 2 [2(D+ 7t? )]3/2 + • • • 



where . . . denote the terms analytic in w. Only the "shear-shear" contribution gives rise 
to the non- analytic |cu| 1//2 term. Using the Kubo formula (I2.23al) . the 0(A) term can be 
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interpreted as the contribution of shear waves to d.c. conductivity, while the |cu| 1/ ' 2 term is 
the cutoff- independent contribution of shear waves to a(u>). The cutoff A is the momentum 
scale at which hydrodynamics ceases being a good effective description. Its value can be 
estimated from the condition that one-derivative terms in the hydrodynamic expansion are 
comparable to zero-derivative terms, A"f v ~ 1. The non-analytic \uj\ 1 / 2 term corresponds to 
the i -3 / 2 fall-off of the correlation function in real time. This phenomenon of "long-time tails" 
in correlation functions of conserved currents has been known for a long time [IS], and, as is 
evident from the derivation, is not specific to relativistic hydrodynamics. In d = 2, one finds 
at zero momentum and small frequency 

k=0) = 0(ln A) + S '^g^i In <^±^ + 0(u » H A>,J)) 



'2 8ttw 2 (D+ 1v 



Again, the O(lnA) term can be interpreted as the contribution of shear waves to d.c. conduc- 
tivity, while the ln(A 2 /|u;|) term can be interpreted in the sense of the renormalization group 
as the logarithmic contribution of shear waves to the "running conductivity" . 
For the correlation function of two stress tensors one finds similar expressions, 



T 



2 



Gf T >,k=0) = O(A)-tfgA- 7+ (|) 



3\3/2 



\UJ 



3/2 



1/2 

+ . . . , d 



G% T Ju, k=0) = 0(ln A) + E% In ^ + . . . , d 



Again, using the Kubo formula (I2.23bj) . the cutoff-dependent terms can be interpreted as 



contributions of shear and sound waves to shear viscosity. 

In d—2 spatial dimensions, one can define the "running conductivity" a(fi) = a(cu—fi), 
and "running viscosity" r](fi) = f](u=fi). Demanding that the correlation functions do not 
depend on fi gives the following renormalization group equations for cr(fi) and i](fi) 

da T x 2 . 

VjT = ~q ~, — -' 3 ' 2a 

Ojj, 87T aw + r/x 

& = -—-, (3.2b) 
ofi 167r rj 

where x — (dn/dfj,)^ is the equilibrium charge susceptibility, and w = e+p is the equilibrium 
enthalpy density [the sliding scale \i in Eq. f)3.2p is not to be confused with the chemical 
potential]. We can rewrite these equations in terms of dimensionless quantities g v = r]/s and 
g a = oT/ x- 

J-^r~ = ' (3.3a) 

on 8wc g a + g v 

M = -"-I, (3.3b) 
Ofi Loire g v 

where the constant c is defined by the equilibrium entropy density as s = cT 2 , in other words 
c measures the number of degrees of freedom in the system. Note that the right-hand side of 
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gain) 




Figure 3: The flow diagram for the differential equations (13. 3p . pictured in terms of the 
dimensionless variables g v = rj/s and g a = crT/x, where r\ is the shear viscosity, a is the 
charge conductivity, s is the equilibrium density of entropy, and \ is the equilibrium charge 
susceptibility. The arrows indicate the direction of decreasing /i (towards lower frequency). 
The dashed line is a straight line with a unit slope, indicating the infrared asymptotics as 
/i->0. 

these renormalization-group equations is small in the large- N limit, and the running of a and 
f] disappears as iV— s-oo. The flow diagram for equations (I3.3P is shown in Fig. |3j Both rj and 
a grow in the infrared, as happens quite generally in two-dimensional hydrodynamics [U] . In 
the limit ji — > the solutions asymptote to g v = g a , in other words 

1 = fl. (3.4) 

s X 

This shows that in the extreme low-frequency limit, 77 and a cease being independent transport 
coefficients in d — 2. 

3.2 Breakdown of second-order hydrodynamics 

The non-analytic behavior of correlation functions of conserved currents at small frequency 
(proportional to |o;| 1//2 in d=3 and In \u\ in d=2) has implications for the derivative expansion 
in hydrodynamics. As indicated in section ll.3[ the derivative expansion is constructed by 
expanding the energy-momentum tensor and the current order by order in derivatives of the 
hydrodynamic variables T, fi, and m m . For example, the shear viscosity arises as the coefficient 
of the one-derivative tensor a^ v in the expansion of T^ v . At second order in derivatives, there 
are other contributions to T^ u , for example <j^ u d\u x . References [HI Q2J SSI El] classify 
possible second-order terms in relativistic hydrodynamics. 

The second-order terms will modify the linear-response correlation functions of T^ u and 
evaluated in Sec. [2J Repeating the calculation in Section 12.61 changes Eq. f!2.33bj) to [TT] 



^T xy T X y k=0 ) =p-iuv + V T nu 2 + ■■■ (3.5) 



where r n is the Israel-Stewart stress relaxation time [49] . and the dots denote contributions 
from other second-order transport coefficients. The derivative expansion is based on the 
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assumption that higher-order terms are smaller than the lower-order terms, and indeed this 
is true in Eq. (13. 5p . in the hydrodynamic limit u— ?-0. 

We saw that in d=3 spatial dimensions, hydrodynamic fluctuations give rise to a non- 
analytic to 3 / 2 contribution to G§ xyTxy (cu, k=0). In the hydrodynamic limit, this is smaller 
than the first-order contribution 0(u), but is larger than the second-order contribution 0(u 2 ). 
Such non-analytic terms can not be represented by a local derivative expansion, which means 
that the derivative expansion is 3+1 dimensional hydrodynamics breaks down at second order. 
This problem with the derivative expansion has been known for a long time [23], and is not 
specific to relativistic hydrodynamics. See Ref. [26] for a discussion of this problem in the 
context of the quark-gluon plasma in heavy-ion collisions. 

We emphasize that this problem with the derivative expansion in hydrodynamics is not an 
intrinsic problem with hydrodynamic equations themselves. The equations of hydrodynam- 
ics can be studied as partial differential equations on their own right at second (or higher) 
order in the derivative expansion. The problem is that these partial differential equations do 
not necessarily capture the correct behavior of fluids in the hydrodynamic regime. Viewing 
hydrodynamics merely as a set of partial differential equations misses information about hy- 
drodynamic fluctuations which render these equations physically invalid as u— >0. A complete 
description of fluids should supplement the hydrodynamic equations with a prescription of 
how to deal with fluctuations. This is the subject of the next section. 

4 Hydrodynamics as a field theory 

So far, we looked at two methods of computing correlation functions of conserved densities in 
the hydrodynamic regime: the standard linear-response formulation in Sections 12.41 12.51 and 
the variational formulation in Section 12.61 There is another way to look at equilibrium real- 
time correlation functions, where the dissipation is attributed to small-scale random stresses 
and random currents in thermal equilibrium [501111]. This formalism has an added advantage 
that the interactions of the hydrodynamic modes can be treated in a systematic way [SI], so 
that both the linear response results of Section &\ and the interaction results of Section [3] are 
included within the same framework. See for example [521 E21 E3]- The stochastic equations 
involving random stresses and random currents can be converted to a functional integral, 
which allows one to view hydrodynamics not just as a low-energy effective description, but 
as a low-energy effective field theory. We will start in Sections 14.11 and 14.21 with equilibrium 
correlation functions derived from linearized relativistic hydrodynamics, which can be viewed 
as propagators for the hydrodynamic perturbation theory. Then in Section |4"U| we will discuss 
interactions between the hydrodynamic modes in a simplified model of Ref. |22j . 

4.1 Simple diffusion 

Just like we did in Section |5J it is instructive to look at a simple example of charge diffusion 
before studying the full hydrodynamics. Beyond being an interesting example by itself, it will 
also serve as a basic building block for the future perturbative treatment of fluctuations. 
For the simple diffusive process studied in Section 12.11 we have found for the retarded 
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function G^ n (u,k) and the symmetrized function G nn (u,k) 

R D X k 2 2TD X k 2 
G nn {u,k) = — — , G nn {u,k) = , n ,,„ ■ (4.1, 



iu-Dk 2 ' "" v ' ' u 2 + (Dk 

Here D is the diffusion constant, and x is the static charge susceptibility, x = (dn/dfi)^ =0 . 
The two functions are not independent, but rather are related by the fluctuation-dissipation 
theorem (12.141) . One can ask the following question: given the hydrodynamic correlation 
functions (14.11) of the charge density n(t, x), what is the effective action S^n] that would 
produce them? Let us now show how one can reproduce correlation functions (14. ip from a 
functional integral formalism, where the charge density n(t, x) plays the role of a "quantum 
field". 



Stochastic model for diffusion 

The diffusion equation (12.11) is a consequence of the current conservation <9 M J M = 0, where the 
charge density is J° = n, and the constitutive relation for the spatial current is = —Dd^n. 
The diffusion equation is only an effective description at long distance and time scales. In 
order to model the microscopic effects, we introduce random currents in equilibrium, so 
that the linearized constitutive relation for the current takes the form = —Ddi-n + . We 
take the microscopic currents to only have Gaussian short-distance correlations, so that 

(r;(x, t) r fc (x', *')) = C5 ik 5(x-x.')5(t-t') . (4.2) 

where C = 2TDx determines the strength of the "noise" . The diffusion equation becomes a 
Langevin-type equation: 

-d t n + DV 2 n + 9(x,t) = 0, (4.3) 
with a random Gaussian field 9 = —dkVk, 

(0(x, t) 0(x', 0) = -2T X DV 2 5(x-x')5(t-t') . (4.4) 

The angular brackets in (I4.4p denote the average over the noise, in other words 

(...) = Jv9e~ w ^.... 

For the chosen noise correlations given by Eq. (I4.4p . the weight factor is 

W n [9] = ^Jdtd d xd d x' 6(x,t)D(x,x.')6(-x',t) , 

where -D(x, x') satisfies 

-CV 2 J D(x,x / ) = 5(x-x'). 



Effective action for diffusion 

One way to compute the correlation functions is to solve Eq. (14.31) and express the density in 
terms of the noise n = n e (x, t), then take the product of the solutions ng(x.i, ti)n e {y.2, i^)---, 
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and average it over the noise. Alternatively, we can pick out the solutions to Eq. (14. 3[) by 
integrating over all possible n(x, t), with a functional delta function which enforces Eq. (14. 3ft : 



n(xi,ti)n(x 2 ,t 2 )---) = Vn(8(e.o.m.)) Jn(xi,ti) n(x 2 ,i 2 ) 



where "e.o.m." stands for the left-hand side of the "equation of motion" ( 14. 3ft . and J = 
det(5(e.o.m.)/5n) is the Jacobian. For a linear equation such as the diffusion equation, the 
Jacobian does not depend on the fields, and can be absorbed into the integration measure. 
We represent the functional delta function by introducing an auxiliary field p(x, t), namely 
<5(e.o.m.) = J Vp exp (i J (e.o.m.)pj . Now the integral over the noise 6*(x, t) can be performed 
because it is Gaussian, giving the following representation of the correlation functions: 



(n(x 1 ,t 1 )n(x 2 ,t 2 )— ) = VnVpe 



-S eS [p,n] 



n(x 1 ,t 1 )n(x 2 ,t 2 



where the action is 



S e s \Pi n ) — Jdt d d x 



C 



{d iP f - ip (d t n - DV 2 n) 



(4.5) 



This gives the functional integral representation of correlation functions for diffusion. We can 
also integrate out the auxiliary field p to end up with the action which depends on n(x, t) 
only: 



(n(xi,£i) n(x 2 ,t 2 )...} 



Vn e 



-S e a[n] 



n(x 1 ,t 1 ) n(x 2 ,t 2 ) 



The effective action for the density fluctuations is 

S cS [n] = ^Jdtd d xd d x' £(t,x)D(x,x')£(t,x') 



(4.6) 



where B(t,x) = (d t n — DV 2 ti) . This effective action does not involve any auxiliary fields, 
it is both real and time- reversal invariant, however it is non-local in space. One can use 
either one of the actions (14. 5B or (14.61) to compute correlation functions. We emphasize 
that the effective action ( 14.61) describing diffusion of a conserved density is not meant to be 
the classical action from which one obtains the "equation of motion" (12. ip via the standard 
Euler-Lagrange variational method. Rather, it is the effective action that is to be used in 
the generating functional for real-time correlation functions of the conserved charge density 
n(t, x) in thermal equilibrium. 



Correlation functions 

What are the correlation functions that we compute in our effective theory (14. 5ft ? In other 
words, what do they correspond to in the operator formalism? It is natural to identify 
(n(xi, ti) n(x 2 , t 2 )) evaluated in the effective theory ( 14.51) with the symmetrized (or unordered) 
density-density correlation function: 

G nri (ti-t 2 ,xi-x 2 ) = (n(xi,ti)n(x 2 ,t 2 )) . 
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This is consistent with the fact that in the path integral formulation, the order of the fields 
does not matter, in other words (n(xj, fa) n(x 2 , fa)) = (^( x 2, fa) w(xi, fa)). With this iden- 
tification, it also important that the effective action ( 14. 6 p is time- reversal invariant. This is 
because time-reversal invariant microscopic dynamics must give rise to time-reversal invariant 
unordered density-density function. However, the action ( 14.51) need not be time-reversal in- 
variant because it also contains information about the mixed (np) correlation functions, and 
those are not time reversal invariant. 

Finally, note that n(t, x) is a real field, and its correlation functions come from the real 
action, Eq. (I4.6p . Therefore, G nn (t, x) is real, or G nn (uj, k)* = G nn (— u, — k). Time-reversal in- 
variance of the action (I4.6P combined with rotation invariance gives G nn (u, k) = G nn (— u, — k), 
in other words G nn (u, k) must be real. 

Let us evaluate the correlation functions in the effective theory (14. 5p . For a quadratic 
action with several real fields a (x, t), we can Fourier transform the fields and write the 
action as 

s = \f C,k^,k)< k . 

The two-point correlation functions are then given by 



G ab (u,k) = (P 



lab 



We apply this to the effective action (14. 5p and read off the correlation functions: 

2T X Dk 2 



G nn (u,k) 
G np (u, k) 



co 2 + (Dk 2 ) 2 
-1 

oo + iDk 2 ' 



G pn (u,k) = 

u — iDk 1 

G pp (u,k) = 0. 

The G nn correlation function has precisely the expected form for the linear response density- 
density correlation function (14. ip . The mixed correlation functions G np (t,x) = {n(t, x)p(O)) 
and G pn (t,x) = (p(t, x)n(O)) look almost like the retarded and advanced functions, and are 
not independent from G nn . To see this, we perform a variable change in the functional integral 
as p{t, x) = tp(t, x) + % J d d x' -D(x, x')5(x', t), and one finds 

(n(ti,x 1 )p(t 2 ,x 2 )) = i Jd d x' L>(x 2 ,x , )(n(t 1 ,x 1 ) J B(t 2 ,x')) . 
The Fourier transform of this relation is 

iuj + Dk 2 

G np (uj, k) = i G nn (u, k) . (4.7) 

Thus we identify — iDxk 2 G np (oj, k) with G^ n (co,k), and — iDxk 2 G pn (u, k) with G^ n (cu,k). 
The relation (14.71) now becomes the fluctuation-dissipation theorem (I2.14p . 
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4.2 Linearized hydrodynamics 

Now that we have worked out correlation functions for the linear diffusion process, we can 
repeat the calculations for linearized hydrodynamics. For simplicity, let us look at relativistic 
hydrodynamics at zero chemical potential, studied in Section E31 



Stochastic model for linearized hydrodynamics 

The hydrodynamic equations are only an effective description at long distance and time scales. 
In order to model the microscopic effects, we introduce random stresses in equilibrium, so 
that the linearized constitutive relation for the stress tensor takes the form 

Tij = Sij (p + v 2 s 5e) - -y v {d^j + dm - ^ij^kj ~ Jc S v 9 k^k + • 

As before, 7^ = T]/w, 7^ = (/w, and w = e+p is the equilibrium enthalpy density, which 
plays the role of the static susceptibility for momentum density fluctuations. We take the 
microscopic stresses to have Gaussian short-distance correlations, 

(T ij {x,t)T kl (x',t')) = 2Tw 

The hydrodynamic equations become 

d t e = -d k n k , (4.8) 
dm = -v%e + Mij-Kj + x) , (4.9) 

where £j = —djT^ is the momentum density noise, M.^ = 7 J) (V 2 <5j J - — didj) + ■jgdidj, and 
7 S = 7^ + ^j^-'j'n as before. Note that there are no noise terms in the equation of energy 
conservation. The noise £«(£, x) has local Gaussian correlations, 

(&(t,x)^',x')> = CV(x-x')5(t-0 , 
where CV,- = —2TwM.ij. The average over the noise & is performed with the weight 

W*[£] = \ Jdtd d xd d x' U^^D^x')^*'^), 

where (x, x') satisfies 



1ri[ S ik 5ji + 8 u 8 jk - jSijhi ) + Ic^ijhi 



8(*-x!)8(t-t') 



Effective action for linearized hydrodynamics 

As before, we can represent correlation functions with a functional integral over the hydro 
variables, which in this case are e and 7r\ Again, we start with 

(e(ti, x x )... } = yp e I?7r(5(e.o.m.))e(t 1 ,x 1 )... 

where the delta function enforces the hydro equations (14.81) and (14. 9p . We introduce an 
auxiliary field f3 for energy conservation, and auxiliary field A; for momentum conservation. 
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Performing the integral over the noise, we find the following representation of the correlation 
functions: 

(7r fc (fi,Xi)7r,(* 2 ,X2) . . . ) = jveVnVpVX e ~ s ^^ ir k (t u Xi)tt,(* 2 , x 2 ) . . . 
where the effective action is 

S cS [e, vr, (3, A] = Jdtd d x fyiCijXj - i\ k F k - iP(dte+d k n k )^ . (4.10) 

Here Fi = d^i + v 2 s dit — Aiijirj. Again, we can integrate out the auxiliary field Aj, which 
gives a representation of correlation functions with a non-local action, 

(7r fc (*i,xi)7r,(* 2 ,x 2 ) . . . ) = JveVirVp e -s«[ e .^ ^{h, xi)tt^ 2 , x 2 ) . . . 

where 

S cfr [e,vr,/3] = ~ Jdtd d xd d x F i (t,x)D ij (x,x')F j (t,x') - i Jdtd d x /3(d t e+d k ir k ) . (4.11) 

Again, we identify two-point functions of the hydro variables evaluated with the action (14.111) 
with the unordered (or symmetrized) correlation functions of the corresponding operators in 
the microscopic theory. The fact that the noise strength Cy is proportional to M.^ ensures 
that the effective action (14. lip is time-reversal invariant. 

Correlation functions 

We can evaluate the correlation functions either from the action (I4.10p or from the action 
(14. lip . Let's work them out using the action (14. lip . To do so, we express the action in 



terms of the Fourier components of the fields, and express the fluctuation in energy density as 
5e = kiTTi/u. The function Fi(u, k) can now be expressed in terms of the momentum density 
alone, 

/ kk\ kk- 
Fi(u, k) = SijTTjiu, k) , where Sij = A v ( 5^ - j + A s -^f , 

where A,, = (—ico + 7^k 2 ), and A s = (—iu + iv 2 ^- + 7 s k 2 ). Setting A^ to zero gives the 
dispersion relation of the shear mode, while setting A s to zero gives the dispersion relation 
of the sound mode. The action is 



so the correlation function is 



Gw j {u>,V) = {S*- 1 D- 1 S- 1 ) ij . 



This can be straightforwardly evaluated: the inverse of is 

(o-i\ _ 1 fx kjkj \ 1 kikj 



V k 2 J A s k 2 
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and the correlation function becomes 

kkA 2T 7r? wk 2 kikj 2 ls wTk 2 uo 2 
~W) u 2 + ^k 2 ) 2 + k 2 (w 2 - v 2 k 2 ) 2 + (7 s k 2 w) 2 ' 

in agreement with the linear-response result (I3.1bl) found in Section 12.41 Similarly to the 
relation ( 14. 7p between G np and G nn for diffusion, one can derive a relation between G^^ and 
G- Ki7Tj in linearized hydrodynamics. Changing the functional integration variable from A, to 
<f>i as Xi(t, x) = <pi(t, x) + % J d d x'Dij (x, x.')Fj(t, x') one finds 

x l) Aj-(t 2 , X 2 )) = I J i d x' D jk (Xi, x') (7Ti(ti, Xi)F k (t 2 , x')) . 

After using uG nit = k\G vm together with the above expression for G ninj , this gives 

k 2 / u + 27,7k 2 k 2 w 2 — v 2 k 2 + iu^k 2 
4.3 One-loop corrections 

So far we looked at fluctuations in linearized hydrodynamics. One could extend the analysis 
of Section 14.21 to the full non-linear relativistic hydrodynamics with the constitutive rela- 
tions (j!.16p . however the resulting effective action is messy and not illuminating. In order to 
study the interactions of the hydrodynamic modes, it is instructive to look at a toy model 
which is technically simpler, yet still retains the essential features of the full hydrodynamics. 
We will thus consider a toy model, with the following simplifying assumptions: i) we only 
expand the hydrodynamic equations to quadratic order in fluctuations, keeping only non- 
derivative quadratic terms, and ii) we impose the "incompressibility" condition d^i = 0. The 
latter amounts to ignoring the "high-energy" sound mode which has u ~ k, while keeping the 
"low-energy" shear and diffusive modes which have u ~ k 2 . 



G T] 



k) 



G ff . A . (w,k) 



Stochastic model for hydrodynamics 

We look at hydrodynamic fluctuations in an equilibrium state with /i = 0. Imposing the 
condition d^k = eliminates the energy density, and leaves n and iik as the only hydro 
variables. Keeping only non-derivative quadratic terms, the constitutive relations for Jk and 
Tij become: 

J k = -Dd k n + — + r k , 

w 

7T/7T/ 

= p5ij - %{di7Cj + dj-Ki) + + nj , 

where r k is the current noise, r»j is the stress noise, and j v = r]/w. The pressure p needs 
to be expanded to quadratic order in n and n k as well, but these terms will turn out to be 
unimportant, as we will see shortly. The hydrodynamic equations give the following stochastic 
model: 

d t n = DV 2 n - -7r k d k n + 6(t, x) , (4.12a) 
w 

d t 7ii = -dip + 7^V 2 7Tj - — n k d k 7ii + x) . (4.12b) 
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From the discussion of diffusion in Section I4.1[ we take the short- distance noise as 

(0(t,x)0(i / ,x / )) = -C a V 2 5(x-x')5(t-t') , 
(&(t,x)£ i (*',x')> = -C,^V 2 5(x-x')5(t-t'), 

with the noise strength constants C a = 2TD\ = 2Ta, and C v = 2T'~y v w = 2Trj. The 
convective term Hkdkn couples the diffusive mode and the shear mode, while the convective 
term -K^d^i gives rise to self-interactions of the shear mode. The model (I4.12p was studied 
by Forster, Nelson, and Stephen in [22]. Note that Eq. (I4.12b|) is just the non-relativistic 



incompressible Navier-Stokes equation with a random external "stirring force" £j. For a 
treatment which includes the sound mode in non-relativistic hydrodynamics, see [55] . 

The symmetry properties of the correlation functions in the model (I4.12p are the same as in 
the linearized hydrodynamics. The equations of the model are real, so the correlation functions 
of n(t, x) and n^t, x) are real. The correlation functions are also invariant under space-time 
translations, and transform covariantly under rotations. As the order of the fields does not 
matter when averaging over the noise, we have G nn (t, x) = G nn {— t, — x), and G Ki - Kj (t, x) = 
Gn-Ki (—t, — x). Combined with rotation symmetry, this implies that G nn (cu, k) and G^ i7r . (cu, k) 
are real (for real u and k). 

Effective action for hydrodynamics 

We follow the same procedure as in Sections 14.11 and 14.21 to derive the effective action for 
the stochastic model (14.121) . As before, there is an auxiliary field p{t,x) for charge density 
fluctuations, and an auxiliary field A«(t, x) for momentum density fluctuations, which we 
choose to satisfy di\ = 0, to match the condition d^i = 0. This condition implies that 
the pressure drops out of the effective action because the corresponding term in the effective 
action J A, dip is a total derivative. This gives the following representation for correlation 
functions: 

(n(ti,xi)n(* 2 ,x 2 )...) = JvnVirVpVX \j\ e -sMn,^ P ,x] n^, Xl )n(t 2 , x 2 ) . . . , 



where 

Ses[n,7r,p, X] 



J dtd d x ^_^ pV 2 p - ^A,V 2 A, - ipB - iXiF^j . (4.13) 

Here B = (d t n — D'V 2 n + Tc^dkn/w), and Fj = (<9 t 7Tj — 7 r) V 2 7Tj + nkdhTVi/w). The conditions 
diiii = and <9jAj = are implied. Again, one could integrate out the auxiliary fields p and 
Aj in Eq. ( I4.13P to arrive at a real (though non-local) effective action. 

The effective action (14.131) is not real. Rather, it is invariant under complex conjugation 
combined with p — > —p and Aj — > — Aj. This symmetry implies that the "mixed" correlation 
functions G np (t, x), G pn (t, x), G Wi \ (t, x), and Ga^ ^x) are purely imaginary. Combined 
with translation invariance, this implies 

G np (u,k) = -G pn (u,ky, (4.14) 
G^ A >,k) = -G A ^(u;,kr. (4.15) 
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Again, we can show that G np is related to G nn , by repeating the same argument which leads 
to Eq. (14.71) . and the extra term (mr k d k n) vanishes by parity. A similar argument gives 
the relation between G n \ and G^, because the extra term {iti'Kkdk'Kj) vanishes by parity. 
Remembering that G nn (u),]<.) and G ninj (u,]t) are real, we have 

G nn (co,k) = 2T X lraG np (co,k) , (4.16) 
G 7Tl7Tj (u, k) = 2Tw Im G niXj (w, k) . (4.17) 

This means that in order to evaluate the correlation functions G nn and G ni7Tj , we can evaluate 
G np and G mi \ , and take their imaginary parts. It will be simpler to do exactly that in practice. 

The above functional integral representation of correlation functions also contains the 
Jacobian J which arises from the variable change to the fields n and 7r. In general, if there are 
several fields (p a satisfying equations E a [ip] = 0, then the Jacobian is J = det(5E a /5(pb). Note 
that in our model the noise terms only appear additively in the stochastic equations, therefore 
the Jacobian does not depend on the noise, and the noise average can be easily performed 
leading to ( 14.131) . We take our fields as ip a = {n, 7Tj}, and we take as their equations of motion 
E a [ip] the corresponding stochastic equations for n and m. Note that the Jacobian is real 
because both the fields (p a and their equations of motion E a [(p] are real. We can exponentiate 
the Jacobian by introducing the anticommuting ghost fields if) a , ip a as was done by Parisi and 
Sourlas in [56] . 

J = Jvi>Vi)e~ s \ S g = Jdtd d xij a ^^-^ b . 

In our model we will have both the density ghosts ip n , ip n , and the momentum ghosts ipi, ipi. 
For the momentum ghosts, we impose the "incompressibility" condition d^i = 0, = 0. 
We can now include the ghosts in the effective action ( 14.131) . and write it as 



S cS = J dtd d x (£ (2) + C {int) ) , (4.18) 



where is quadratic in the fields, and £( m *) contains the interactions, 

£ (2) = -^pV 2 p-^A,V 2 A 4 -^(^- J DV 2 n)-^(9 t vr 4 -7,V 2 7r J ) 

+ Ud t - 7 ,V 2 )^ + Md t - DV 2 )iP n , (4.19) 



£M) _ _Ap 7r Q n _ L\ n Q n 
W W 

+ —4>i{.d k TXi)4> k + — 4)i7l k d k 4ji + —$Jdin)$i + — ^nTTfc^fcV'n • (4.20) 

w www 

The conditions d^i = 0, di\ = 0, ditpi = 0, = can be imposed by additional Lagrange 
multipliers. These extra auxiliary fields only appear in £( 2 \ but not in £( mi ) ; hence their 
only effect is to modify the propagators, but not the vertices. The action above is suitable 
for calculating correlation functions in the hydrodynamic perturbation theory. 
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The action (I4.18P is invariant under space-time translations, spatial rotations, parity, and 
charge conjugation. The action ( I4.18P is not invariant under time-reversal, nor should it be. 
The action is also invariant under Galilean symmetry under which 



1 1 



7Ti(t, x) — > 7T.j(t, X— Vt) + WU 

n(t, x) — y n(t, x— vt) , 

and all other fields transforming the same way as n(t, x). Further, the action (14. 18j) is invariant 
under the symmetry generated by 

5n = £ip n , 5tp Q = i£p, 

where £ is an infinitesimal anticommuting parameter. This symmetry is quite generic in 
stochastic systems J56j[57]. It is an analogue of the BRST symmetry, and implies the relations 



^n^n(^' X ) ~~ —iG rip (t,X.) , 

GWj (*> x ) = -* G-KiXj (t, x) • 

For the diagram technique, we read off the propagators from C^ 2 \ and the vertices from £( m *) 
The propagators were already calculated in Sections 14.11 - 14.21 and we summarize them here 
associating a line with each one: 



C k 2 



1 



2^2 



Gn P (u, k) — ^ ijj^ 



1 

u - iDk 2 



Gppi^i k) — 'WWWWWWWW 



I J 



i 3 
nrtsirsisisirsirsirsirsisis^ 



G° AiA >,k) = 

The vertices are shown in Figure HI There are also ghost propagators and ghost vertices which 
we didn't write down, but they can be easily read off from the action (I4.18P . 

Correction to the momentum density propagator 

At one loop, there are five connected diagrams potentially contributing to G. Ki \ j , three of which 
are shown in Figure [51 and the remaining two contain ghost loops. The first two diagrams 
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Figure 4: The vertices for the effective action (j4.13j) . There is an overall sign depending on 
whether the momentum is flowing in or out of the vertex. 




Figure 5: Three connected one-loop diagrams potentially contributing to G n \. In addition 
to the shown diagrams, there are two more diagrams with ghost loops which are completely 
analogous to the diagrams in the first row. 

shown in the figure vanish. To see this, let u and k be the external frequency and momentum, 
and z and q the frequency and momentum running in the loop. The first diagram will have 
G np and G pn in the loop, and is proportional to 

r dz -l i 

J 2^ z + iDq 2 (u-z) - iD(k-q) 2 ' 

In the plane of complex z, both poles are in the lower half plane, and as a result the integral 
is zero. Exactly the same argument applies to the second diagram in Figure and to the 
diagrams with ghost loops. Thus it is only the last diagram in Figure which contributes. 
Let us call the amputated part of the last diagram E mn (u;,k), i.e. we write the connected 
one-loop contribution in momentum space as 

By rotation invariance, E mn must have the form E mn = 5 mn Ei + A; m fc n E 2 . Because Gj. A . is 
transverse, it is only Ex that contributes. By summing all such one-loop diagrams, we find 
the correlation function as G = G°(l — ExG ) -1 , which gives 

If Ei(w, k) is proportional to k 2 for small external momenta, its imaginary part can be in- 
terpreted as a correction to 7^, or equivalently as a correction to viscosity. The correlation 
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GW>, k^O) = ( 8 i3 - ^ ) , (4.21) 



function of momentum density then becomes 

2Twj v {u)k 2 f x 

L0 2 +(j v {u) k 2 

where 7,7(0;) = j v + ImE^w, k— >-0)/k 2 . Explicitly, the loop integral is 

v, / , ^ k 2 f d d q ( d 2 -d-2 2 2i lv q 2 \ , t nnS 

where s = w/T is the equilibrium density of entropy. In the large- N limit, both the viscosity 
and the entropy density grow with N, so that 7^ = 0(1), and 1/s can be treated as a small 
parameter, suppressing the fluctuation corrections. In d = 3 spatial dimensions the integral 
(I4.22p is linearly divergent at large momenta. Introducing a large-momentum cutoff, we find 



Re£i(o;,k^O) = v ' 1 ° ; ; , ImSi^k-^O) = C A k 2 



30tts (4 7r ,) 3 / 2 ' iV ' 30tts (4 7?? ) 3 / 2 ' 

The constant Ca is proportional to the cutoff and to 1/s, but is u and k independent. The 
imaginary part of £1(0;, k) can be interpreted as a correction to 7^, producing a frequency- 
dependent viscous damping coefficient 7,7(0;) = 7,7 + Ca — 3§^\/M (47»?) _3 ^ 2 - We define the 
"renormalized" viscous damping constant as 

7„ = % + <-' A- 

We need to express the physically measurable 7,7(0;) in terms of finite 7^ en . Following the 
standard renormalization procedure, to this order in 1/s we replace j v with 7^ en , which gives 



23 i\uj\ 

The renormalized viscous damping constant in d = 3 is given by the usual Kubo for- 
mula (12.22bj) . as u — > (recall that rj = 7,7 w). For example, taking the spatial momentum 
along z, we have 

l7 = lim lim ^G Wm (u, k) . (4.24) 
' 21 w w->ok->o k z 

In d=2 spatial dimensions, the integral (I4.22p is logarithmically divergent at large mo- 
menta. Introducing a large-momentum cutoff, we find 



Re E x (o;, k— )-0) = — — , ImEi(u;, k-)>0) = — . 

32s 27^ 1d7ts 27^ 

Again, the imaginary part of £1(0;, k) can be interpreted as a correction to j v , producing a 
frequency- dependent viscous damping constant, 7,7(0;) = + ^C^/M)- ^ we interpret 

the logarithmic divergence in the renormalizat ion-group sense, we can define the "renormal- 
ized" damping constant 7™ 11 as 7,7(0;) evaluated at some frequency scale uj = /i, 

1 1 , A 

^ = ^ + mTs2^ l %- 
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AA/WV 



Figure 6: The connected one- loop diagram contributing to G 



np- 



To this order in 1/s we replace 1/7^ in the right-hand side with l/7^ en , and express the "bare" 
parameter 7^ in terms of 7^ en . The "bare" parameter j v knows nothing about the arbitrary 
frequency scale \x at which we choose to define the viscosity. Hence d^/d/j, = 0, which implies 

dY en 1 1 

a/x 327rs 7^ en 

In other words, in order for the frequency-dependent viscosity coefficient 7^(w) in the corre- 
lation function (I4.2ip to be independent of the arbitrary scale /j,, the renormalized viscosity 
must depend on the arbitrary frequency scale //, at which we choose to measure it, so that 
7^ en (/x) grows in the infrared, proportional to A/ln(l/yu). The viscous damping coefficient 
-f v = r]/w is 

7 ( w ) = + — In ^ , 

and the correlation function of momentum density ( 14.2ip does not depend on fi. The correction 
to differs from the one computed in Section I37L1 because the present "incompressible" model 
ignores sound waves. 

Correction to the charge density propagator 

Let us now calculate the correction to the charge density correlation function. At one loop, 
there is only one connected diagram that contributes to G np , as shown in Figure El Let us 
call the amputated part of the G np diagram k). By summing all such diagrams, we find 
the corrected correlation function: 

oj + iDk 1 + h(oj, k) 

For small external momenta, E(tu,k) is proportional to k 2 , and we obtain the correlation 
function of charge density as 

CU«,k->0)= ZTXDM* (4.27) 
co 2 + \ D(u) k 2 J 

where D(u) = D + ImS(w, k— >0)/k. 2 . Explicitly, S(o;, k) is given by 

S(w,k) = ^ I ^G^.(z,q)G;( U -,,k-q)ft- gi )fc r (4.28) 
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In the expression (14.281) . the g» factor does not contribute because G° j7r . is transverse. The 
frequency integral can be evaluated by closing the contour in the lower half-plane. The 
remaining momentum integral simplifies when k — > 0, and one finds 

k^O) = -*±± [** / , 2 . (4.29) 

V ; s d J (2vr) d to + i{l v +D)c? V ; 

In d = 3 spatial dimensions, the integral (I4.29P is linearly divergent at large momenta. Intro- 
ducing a large-momentum cutoff, we find 



Re£(c,k^0) = -^ r v o , ImS(c,k->0) = C A k 2 k " X M 



3vrs [2( 7j? +L>)] 3 / 2 ' v ' 3tts [2( 7 ,+D)] 3 / 2 ■ 

The constant Ca is proportional to the cutoff and to 1/s, but is u and k independent. The 
imaginary part of E(a>, k) can be interpreted as a correction to D, producing a frequency- 
dependent diffusion coefficient D(u) = D + Ca — t^a/M [2( 7j? +-D)]~ 3 / 2 . We define the 
"renormalized" diffusion constant as 

D Icn = D + C A . 

We need to express the physically measurable D{uS) in terms of finite D ren . Following the 
standard renormalization procedure, to this order in 1/s we replace ( / y v +D) in the right-hand 
side with (Y v en +D ven ), which gives 



D(u) = D rcn - — — v H N1 , /0 . (4.30) 

v ; 3vrs [2(7^ n +D ren )] 3 / 2 v ; 

The renormalized diffusion constant in d = 3 is given by the usual Kubo formula (I2.22a[) . as 
co — > (recall that a = Dx), 

1 u? 

D rcn = Um Um q / jx / 4 31 x 

2Tx^ok^ok 2 v ; y ' 

In d = 2 spatial dimensions, the integral (I4.29P is logarithmically divergent at large mo- 
menta. Introducing a large-momentum cutoff, we find 

. . k 2 sign(w) . . k 2 ln(A/M) 

Re E(u, k— >0) = -— , ImE W ,M = « ^TT^ • 

16s 'Jn+D Sirs ^+1) 

The imaginary part of k) can be interpreted as a correction to .D, producing a frequency- 
dependent diffusion coefficient, D(uS) = -D+g^ 7 ^ D ln(A/|a;|). If we interpret the logarithmic 
divergence in the renormalization-group sense, we can define the "renormalized" diffusion 
constant D Ten as D(u) evaluated at some frequency scale uj = /i, 

D ron = D + — — In - . 

87rs 7r) -|--D n 

To this order in 1/s we can replace {p/ n +D) in the right-hand side with (7^ en +D ren ), and 
express the "bare" diffusion constant D in terms of D ien and 7 ^ en . Note that the "bare" 
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D(/i)V7 




Figure 7: The flow diagram for the differential equations (I4.33p . pictured in terms of the 
dimensionless variables 7] / {Ty/s) and D y/s, where i] is the shear viscosity, D is the diffusion 
constant, and s is the equilibrium density of entropy. The arrows indicate the direction 
of decreasing fi (towards lower frequency). The dashed line is a straight line with a slope 
(y/Vf— 1)/2, indicating the asymptotics as /x — > 0. 



parameter D knows nothing about the arbitrary frequency scale /j at which we choose to 
define the diffusion constant. Hence dD/dfi = 0, which implies 

dD ven 11 . . 

A*-* — = — — . 4.32 

dfi 8-ks 7^ en + J D ren v ' 

In other words, in order for the frequency-dependent diffusion coefficient D(u>) in the corre- 
lation function (I4.27P to be independent of the arbitrary scale the renormalized diffusion 
constant must depend on the arbitrary frequency scale /i, at which we choose to measure it. 
The diffusion coefficient is 

D(u) = D rcn (n) + — , N 1 - ~. t In ^ , 

and the correlation function of charge density (I4.27P does not depend on [i. 



Scale dependence of the diffusion constant and viscosity in d = 2 

We have found that the viscosity and the diffusion constant are running, scale-dependent 
parameters in d = 2 spatial dimensions, satisfying renormalization-group equations (I4.25P 
and (I4T32D . 

frC 1 L dDTCn - 1 1 (4 33) 

where D is the charge diffusion constant, 7^ = 77/ (e + p) is the shear mode damping con- 
stant, and fj, is the frequency at which we choose to measure the transport coefficients at 
zero momentum. These equations are almost identical to the equations (13. 2p derived earlier, 
except for a factor of 2 in the equation for 7^. This is because the present model only keeps 
shear waves, but neglects sound waves, which contributed to the renormalization of the shear 
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viscosity in Section 13.11 In the low-frequency limit \i — V 0, Eq. (I4.33P give the following scale 
dependence of the transport coefficients: 



The flow diagram for Eq. (14. 33p is shown in Figure [7J One can see that in the extreme 
low-frequency limit, the renormalized shear viscosity and the diffusion constant (or charge 
conductivity) cease being independent transport parameters, but are rather related by 

DT= Vi7-iv ^ 156 v_ (434) 



5 Discussion 



These lectures mainly focused on response functions of conserved densities in relativistic 
fluids. In linear response theory, hydrodynamic correlation functions may be evaluated by 
several methods: the canonical approach of Sec. I2.2[ the variational approach of Sec. 12. 6[ 
or the stochastic approach of Sec. 14.21 In non-relativistic fluids, the corresponding response 
functions can be measured using light scattering [58]; relativistic fluids are far less common. 
It would be wonderful to see relativistic hydrodynamics established experimentally as the 
effective description for many-body quantum systems with emergent Lorentz symmetry, such 
as those discussed in Ref. [2J. 

There are a number of subjects related to hydrodynamic fluctuations in relativistic fluids 
that we have not touched on. These include hydrodynamic three- and four-point functions 
[59| loP] , hydrodynamics of anomalous currents [151 ED E2J E3] , hydrodynamics of relativistic 
superfluids [HH H7J E5J, hydrodynamics of theories with broken parity [THl H7J, relativistic 
magneto- hydrodynamics (55] , and the relation of the noise to the closed time path formalism 
in quantum field theory [57] . 

As we have seen in Section 13. 2\ fluctuation effects render second-order hydrodynamics 
invalid. On the other hand, second-order hydrodynamics was introduced in the first place 
in order to eliminate the acausality of the first-order hydrodynamics associated with short- 
wavelength modes. A complete hydrodynamic theory must contain both the derivative expan- 
sion, and a systematic procedure to treat the fluctuations. 

In these lectures, we tried to emphasize that hydrodynamics should be viewed as more than 
a study of the classical hydrodynamic equations. While many questions in hydrodynamics 
can be answered by viewing it as a classical field theory, there are aspects of hydrodynamics 
which are more quantum field-theoretic in nature. The derivative expansion of the constitutive 
relations in hydrodynamics is not unike the derivative expansion in effective field theory, the 
hydrodynamic equations are analogous to the effective Lagrangian, and the effects of the 
hydrodynamic fluctuations are analogous to quantum loop corrections. A notable difference 
between hydrodynamics and effective field theory is that at the moment there is no simple 
and systematic procedure to reformulate hydrodynamics as a field theory, where the derivative 
expansion can be applied directly to the effective action^ A naive attempt at deriving the 

11 Formulating an effective action for hydrodynamics appears an easier problem if one ignores dissipation: 
see for example Refs. [H |69l EQl ED HH Eg] . 
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effective action for relativistic hydrodynamics along the lines of Section H] leads to an action 
which is not particularly simple or illumminating. We leave the problem of the effective action 
for dissipative relativistic hydrodynamics including the questions of frame invariance and the 
derivative expansion for future work. 
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